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ABSTRACT 


Accurate simulation of helicopters in an at-sea-shipboard environment is desired to 
provide realistic operating envelopes without incurring the enormous cost of real-time 
flight tests. This study examines the simulation of rotorcraft in turbulent flow by look- 
ing at previous attempts at helicopter-ship interfacing, current efforts in this area, and 
What will be needed in the future. Part of this study is devoted to the construction of 
an analytic model of the “tunnel strike” problem of the CH-46 Sea Knight helicopter that 
is based on measurements made over the flight deck of a model ship. A computer model 
was constructed, with the aim of modeling the “tunnel strike” during engagement aboard 
AOR type ships. The remainder of the study is concerned with the simulation of the 
motion of a helicopter in the turbulent wake of a DD-963 class ship. Results show that 
a sixth order transfer function can filter white noise to accurately model the turbulence 
spectra at specific points along a helicopter glide path in the wake. While a tunnel strike 
could not be successfully modeled using DYSCO software, a simple blade-element pro- 
gram was developed to show the aerodynamic forces on the rotor blades in a specific 
flow field over an AOR class ship flight deck. That program shows the location and 
magnitude of the aerodynamic forces contributing to the flapping of the rotor blades, 


which results in the rotor blades impacting the fuselage. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 


additional verification is at the risk of the user. 
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I. INTRODUCTION 


Interface flight tests are used to determine the operating envelopes for various heli- 
copters on the various classes of U.S. Navy ships. The high cost of ship-helicopter 
interface flight tests indicate the need for a faster and less expensive method to determine 
the safe operating envelopes of the fleet’s helicopters on aviation capable ships. One 
solution to this is to attempt to develop an accurate, real-time computer simulation of 
various ship-helicopter combinations [Ref. 1]. In order to accomplish this a computer 
model must include a sufficiently accurate math model of the dynamics and aerodyna- 
mics of the helicopter, and the operating environment. Several attempts have been made 
to model the aerodynamic wake turbulence for various class ships, all of which have re- 
sulted in simulations that bear no relation to the real flight. One such attempt, made 
by Fortenbaugh, was the use of Strouhal number scaling to extrapolate the ship air-wake 
model of a FF-1052 for a DD-963 class [Ref. 2]. 

For aircraft simulations in land based flight, the current MILSPEC dictates the use 
of a Dryden stochastic model [Ref. 3]. Since the wake of each class of ship 1s different 
no such standardization is possible. This turbulence has a great influence on the oper- 
ation of helicopters aboard Navy ships, and needs to be modeled correctly in order to 
provide an accurate simulation for pilot training. The greatest concerns at the present 
time are methods used to combine ship air-wake turbulence models, such as those de- 
veloped by Fortenbaugh [Ref. 2] and Hanson [Ref. 4] with an accurate helicopter math 
model that gives a true representation in real-time of the helicopter’s motion in all flight 
profiles. This is extremely difficult to do in any quantitative way, considering the nu- 
merous stability derivatives present in a six degree-of-freedom helicopter model. 

Another area of particular interest at the present time is the rotor engage/disengage 
envelope of Navy/Marine helicopters on aviation capable ships at sea. This interest is 
due to the occurrence of rotor blade strikes to the fuselage during rotor engagement and 
disengagement. This problem has occurred on the AOR and LHA class ships frequently 
enough to warrant investigations into their cause. In order to complement completed 
studies of airwake aerodynamics of the flight deck of the AOR class ships [Ref. 5], this 
study will concentrate on the AOR class ship. Current limitations are illustrated in 
Figure 1 on page 2. 


AOR CLASS SHIPS 


ROTOR ENGAGE/DISENGAGE 


WIND LIMITS 


All Wind Azlmuths Relative 
to Shlip’s Centerline/ 


Alrcraft Parallel 
To Ship's 
Centerline 


345 to 015/35 
016 to 010/30 
041 to 180/45 
181 to 235/25 
236 to 320/15 
321 to 344/25 


Alrcraft Parallel 


To Port-To-Star- 


board Landing 
Lineup Line 


340 to 005/45 
006 to 035/35 
036 to 050/30 
051 to 070/20 
071 to 080/30 
081 to 180/45 
181 to 235/25 


MaxIlmum Wind 
Speed (Knot) 


Alrcraft Parallel 
To Starboard- 
To-Port Landing 
Lineup Line 


345 to 005/45 
006 to 025/40 
026 to 040/30 
041 to 180/45 
181 to 255/20 
256 to 325/25 
326 to 344/40 


236 to 310/15 
311 to 339/35 





Figure 1. CH-46 Rotor Engage/Disengage Limits for AOR Class Ship. 
Source: NAVAIR AI-H46D-NFM-000, H-46D NATOPS Manual 


Due to both the atmospheric turbulence and that caused by the airflow around the var- 
ious superstructure elements of the ship, a highly turbulent recirculating flow exists on 
the helicopter landing deck, typically located on the aft end of the ship. Previous work 
by Anderson has attempted to measure these flows in a wind tunnel environment that 
simulates the air-wake turbulence of a DD-963 class ship [Ref. 6]. These unsteady tur- 
bulent flows restrict the starting and stopping of rotor systems due to excessive flapping 
of the blades. One of the most severe restrictions apply to the CH-46 Sea Knight tan- 
dem rotor helicopter, built by Boeing Helicopter Company (Boeing Vertol). The “tunnel 
strike” occurs when the blades of the aft rotor head bend down far enough to strike the 
fuselage and housing or “tunnel” that covers the interconnecting shaft between the aft 
and forward rotor systems, near the ‘Y= 180° position, with ‘’=0° located along the 


longitudinal axis as shown in Figure 2 on page 3. 





Figure 2. CH-46 Sea Knight Relative Blade Position. 

Source: NAVAIR Al-H46D-NFM-000, H-46D NATOPS Manual 

Analytical studies by Leone at Boeing have indicated a possible source of the problem, 
but no solution or safe operating envelope has been found without using actual flight 
tests [Ref. 7]. If a computer model of the “tunnel strike” could be made using either 
math models or physical models of various ships, determination of safe envelopes could 
be done faster and cheaper. One possible solution is to construct a computer model of 
the H-46 and expose it to a physical model of the flow over the flight deck in an attempt 
to create tunnel strikes constructed from current mishap data. One potential modeling 
tool is DYSCO, a dynamic system coupling program produced by Kaman Aerospace. 
The details of this system will be discussed later. By introducing velocities obtained from 
Computational Fluid Dynamics (CFD) programs or wind tunnel tests, it may be possible 
to simulate strikes for various helicopter positions on the flight deck, as well as for dif- 
ferent class ships. A limitation of using the velocities from CFD programs is that the 
power spectra of the velocity fluctuations cannot be predicted at present, and thus can- 
not account for other than mean airflow velocities. If successful, such a method could 
replace the expensive and time consunung flight tests now used to determine safe rotor 


engagement/disengagement envelopes. 


Thus, the goals of this project are twofold. First, to investigate the steps that would 
need to be taken to incorporate a stochastic ship air-wake turbulence model into a heli- 
copter math model capable of reacting accurately to the multi-directional, multi-valued 
velocities and produce a real-time simulation with sufficiently fast computation time. 
Second, to use the wind component velocities obtained from a CFD program (not yet 
developed) or from tunnel tests at the Naval Postgraduate School in an analytic code to 
duplicate actual tunnel strike incidences, and then investigate with this model possible 


safe engagement/disengagement envelopes for specific ship classes. 


A. PREVIOUS STUDIES 
1. Air Wake Studies 

Studies of the airwakes of ships were done at Boeing Vertol Company in 1976 
for a FF 1052 class frigate and 1980 for a DD-963 Class destroyer. These air-wake 
models were obtained using scaled models of the ships in a uniform-flow wind tunnel 
[Ref. 8} [Ref: 9}. The details of how the models were obtained are explained by 
Anderson [Ref. 6: p. 4.]. The FF-1052 model was used in 1978 by Fortenbaugh, in an 
attempt to generate data for the DD 963 class destroyer and used Strouhal scaling to 
relate the frequencies. This scaling procedure is valid for relationships between different 
sizes of the same ship, not for ships of different configurations. 

Based on the recommendation of Fortenbaugh, the airwake of a scale model of 
the DD 963 class destroyer was mapped in a wind tunnel by Boeing in 1980 [Ref. 9]. 
Again, the same assumption of no environmental turbulent boundary layer was used. 

A study by Nave produced several algorithms with which to mathematically 
model the ships’ airwakes for incorporation into a motion simulator [Ref. 10: p. 6.]. 
As with all simulators studied, Nave’s turbulence model assumes a point mass aero- 
dynamic model in which all aerodynamic forces and moments are calculated from the 
mean value of airspeed, angle of attack, and sideslip at the aircraft center of gravity 
[Ref. 10: p. 16 }. This is due to the availability of single point data only, as noted by 
Anderson [Ref. 6: p. 2.], since no multi-point model for helicopters exists. Using a 
second order filter equation, Nave investigated several variations of damping and filter 
frequency. The filter helps shape the straight line Gaussian noise to the values required 
to represent the power spectral density of the velocity fluctuations in the ship’s airwake. 
The form of the filter is: 


K(S + a) 
LS ieee 
S° + 20@,S + w,, 


OY 


59 
frequency, a is the minimum phase zero of the transfer function and V is the airflow 





using the baseline of w, = and ¢ = 0.4, where K 1s the filter gain, w, is the system 
mean velocity, as suggested by Garnett. [Ref. 10: p. 28.] (Ref. 8] This configuration 
matched the characteristics of the measured data. Unfortunately, these relations are 
derived from wind tunnel data in a uniform flow, and without the use of an environ- 
mental boundary layer and they must be viewed with some skepticisin. 

Hanson [Ref. 4] also developed a turbulence model, based initially on the 
airwake data from the DD-963 model developed by Fortenbaugh. Hanson refined 
Nave’s analysis of processing the output of the random number generators prior to being 
input into first order shaping filters. This was done by generating the random numbers 
at intervals other than the time frame of the simulation and linearly interpolating the 
random number sequence for time values between the successive calls to the random 
number generator, in order to give a smoothly varying number sequence. This method 
is shown in Figure 3, where 7, is the interpolated Gaussian white noise input, and vy, is 
the filtered white noise output (Ref. 4: p.11] 


n 
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Figure 3. Block Diagrain of Linear Interpolation Scheme 
Source: Hanson, G.D., Airwake Analysis 


Nave’s method [Ref. 10] produced turbulence levels that were too high, ac- 
cording to pilot reports on the flight characteristics of the model as installed in the 
NASA Ames H-2 simulator model, and thus required modifications. These modifica- 
tions included changing the means and variances of the FF-1052 data base to a DD-963 
data base, elimination of the Strouhal scaling relations with development of a DD-963 
wind tunnel data base [Ref. 9], and changing the extrapolation of the shaping functions 
to exponential versus the cosine type functions used by Nave (Ref. 4: p.31.]. 

Anderson’s study of DD-963 Airwake Turbulence [Ref. 6] at the Naval Post- 
graduate School included measurements along three typical helicopter landing paths, but 
the data was defective due to a hot-wire system malfunction. These data points have 
since been retaken and are discussed in this study and presented in Reference 11. 

2. Helicopter Models 
a. GENHEL 

One helicopter simulator in use today is GENHEL, developed by Sikorsky 
Helicopters, and run on a CDC 7600 machine. This is a six degree of freedom (DOF) 
system, and was used in the development of the UH-60 Black Hawk. The flow diagram 
for this simulator is illustrated in Appendix C. The atmospheric turbulence model used 
is the Dryden model, as required by current military specifications [Ref 3: p. 45.]. This 
model is represented by the algorithm: 


OSU) age 
ee hee (2nn) 
where n=n =, n being the frequency, L turbulence length scale and U the mean com- 
ponent velocity. This function is used frequently because of its simplicity; however, the 
accuracy of this function falls off with increasing frequency [Ref. I: p. 14.].[Ref. 12: 
p.2]. Inputs required for the GENHEL program include: 
Aircraft velocity relative to the airmass 


Longitudinal turbulence length scale 


Root mean square (rms) gust velocities 


These inputs provide the environment for calculating the effects of the turbulence on the 
fuselage and rotor blades. The blade forces are found by strip theory, summing up the 
contributions of each blade for forces and moments. The rotor forces and moments are 


obtained in a wind-hub coordinate system, then transformed to a hub-body system. 


These forces and moments are then given in reference to the body CG by transforming 
all forces and moments to body axes.[Ref. 12: p. 4.] 
b. Parallel Processing 

One of the problems with introducing non-steady flows to rotary wing 
aerodynamics is the large amount of computation necessary to achieve an accurate sol- 
ution. Most of these effects require implementation of the Navier-Stokes equations, 
which are expensive to solve. Such a solver was introduced by Wake and Sankar, but 
full scale Navier-Stokes are not yet feasible on the CRAY-2 supercomputer due to the 
extremely fine grid system required to resolve the boundary layer [Ref. 13]. Since 
Supercomputers are expensive, and the Navier-Stokes equations cannot be solved in real 
time, both are impractical for simulator use. 

One approach is to simplify the entire model. This has been accomplished 
by Messerschmitt-Belkow-Blohm GmbH (MBB) in West Germany. The rotor aero- 
dynamic model is based on blade-element theory, including the effects of compressibility, 
stall, and reverse flow effects. Also included are rotor downwash (by modified momen- 
tum theory) and the influence of ground effect on the rotor. This, along with rotor 
flapping dynamics are considered for each blade separately. All blade forces are com- 
puted simultaneously by a separate processor for each blade, and when each processor 
is finished, the total rotor forces are computed. See Figure 4 on page 8 for the flow di- 
agram of this process. The total simulation model is presented in Figure 5.[Ref. 14] 


Simulation Integration Step At Rowen latlonniode: 


Equivalent System for 
Flapping for the 
act. Rotor Speed 


Starting of the 4 Parallelprocessors 





Rotor Speed 


Biade 1 Blade 2 Blade 3 Blade 4 





; 





Wait for the 4 Parallelprocessors to Finish 


Total Rotor Forces 


Figure 4. Parallel Calculation of Blade Dynamics (4 blades) 
Source:Huber, H., AGARD CP 359 


While modeling the unsteady aerodynamics is difficult, modeling a helicop- 
ter in atmospheric turbulence only complicates the problem, and dims the prospect for 
a real-time solution. MBB, using their parallel processors, introduce discrete gusts and 
continuous stochastic disturbances in accordance with MIL-F-8785C [Ref. 3: p. 45.]. 
The stochastic disturbances are simulated by passing a random signal (white noise) 
through a first order lag function, and adjusting the filter gain and time constant to 
produce a power spectrum match with the Dryden model [Ref. 15: p. 68.]. That lag 


function 1s: 


A 


2 ears 


where A is the filter gain, and B 1s the location of the single pole of the system. Taking 
into consideration the Dryden model for the horizontal gust component, the horizontal 


turbulence power spectra model would be: 


/ 20 ik ! 
FAS) = a Ary 1+ L,,|vs 


where F,(S) is the power spectra, with the sampling frequency At and the factor “a” to 





adjust the standard deviation of the noise source. Here, L, is the horizontal length scale, 
o, is the horizontal variance of the noise, and v is the magnitude of the velocity vector 
of the wind disturbance. This method is employed in this study to obtain suitable 
models for the power spectra in the wake of a ship. [Ref. 15: p. 69.] 


B. A PREVIOUS TUNNEL STRIKE STUDY 

Early investigation of the excessive flapping of CH-46 rotor blades was conducted 
in 1964 by Peter Leone [Ref. 7]. He studied the transient aero-elastic response of the 
helicopter blade, and the excessive flapping and droop-stop impact. This study was done 
with the rotor RPM at one hundred percent, with a step input to the rotor system from 
the cyclic pitch controls. Using the Myklestad finite difference equations, Leone calcu- 
lated the fundamental and first modes in flapping [Ref. 7: p.36.]. This assumes only 
cantilever flapping motion, without consideration of coupling with the lead-lag and 
torsion modes. This approach also used only mean wind flows, and did not account for 
discrete or random wind gusts. In 1982, at the request of NAVAIR, Boeing Helicopter 
conducted an analytical study of excessive flapping at low rotor RPM [Ref. 16]. This 
report noted that upflow through the rotor disk contributed to the high bending de- 


| Pilot Cues | 
Sight/instrument 
iDlsplays/Buffetting 



















Computer Generated 
Image System 













Aerodynamic 


ee 
| | [| | Fuselage 


Pilot Iinteractlo 
Control Input 


Relative Motion 
in the 
Eye Point System 


— 
| | __|Main Rotor 


External Body 
Axis Forces 


ro 


Motlon in the | 
7 | Interface 
Flight 


Suatem ! Computer 
Equations 
Numerical | 

| Integration | Wind 






Avionics 


Engine 





Atmosphere 
4 
a ed 


Figure S. Total Simulation Model (MBB) 
Source: Huber, IL, AGARD CP 359 


flections with the blades. Using a company program labeled as LO2 and a RPM accel- 


eration rate suniliar to that on page 31, the maximum bending deflections to hit the 
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fuselage occurred at 16 RPM, which agree with tunnel strikes reported before and after 
Leone’s study. Leone’s assumptions of the angle of attack of the wind through the rotor 
system were based on wind tunnel flow visualization techniques of an oil rig [Ref. 7: 
p.5.]. Figure 6 on page 12, Figure 7 on page 13, and Figure 8 on page 14 illustrate the 
results of Leone’s investigation. Figure 6 shows the position of the rotor disk with re- 
spect to the flight deck and incoming wind vector, and the rationale for the angle of at- 
tack of 15°. Figure 7 illustrates the upward blade displacement as a function of main 
rotor RPM, while Figure 8 shows the downward blade displacement as a function of 
main rotor RPM. Using this angle of attack, he broke down the wind to vertical and 
horizontal components, addressing only the mean velocity for each component. Using 
the assumption that this vertical flow from the side of the ship would cause the rotor 
blade to “sail” up until it impacted the upper flap restraint, he computed the blade tip 
deflection using the first bending mode. 

Unfortunately, the geometry of the oil rig is a poor representation of the landing 
deck of a ship. Thus, the low (15°) angle of attack assumption made here causes Leone's 
results to be suspect. While there is no argument of the existence of an upflow as the 
wind meets the side of the ship, flow visualizations by Johns [Ref. 17] indicate that 
turbulence levels exist that would require the use of more than just the mean velocities 
to compute the effects of the vertical relative wind component on the exposed rotor 
blade. The behavior of the upflows on the AOR were investigated by Rhoades [Ref. 
5] who made measurements at four points around the blade-tip trajectory of the aft 
CH-46 rotor for a range of ship-wind relative velocity vectors. Mean values of the ve- 
locities are used in Chapter III due to the unavailability of more suitable values and the 


simplicity of the model developed. 
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Il. THEORY OF TUNNEL STRIKE AND AIR-WAKE TURBULENCE 
MODELS 


A. DYNAMIC SYSTEM COUPLER PROGRAM (DYSCO) 

The program DYSCO is used in this study to model the flapping behavior of the 
main rotor blades of the CH-46 helicopter. DYSCO 1s installed on a Digital Electronics 
Corporation VAX 2000. It is an interactive computer program which allows the user 
to model the dynamic and aerodynamic behavior of rotorcraft and other aerospace 
structures [Ref. 18: p.1]. This is done using component modules to represent the 
second-order ordinary differential equations which may have constant, time-dependant, 
or non-linear coefficients. For the CH-46 Tunnel Strike model, the following compo- 
nents were used: 

CFM2 Two Dimensional Moda! Fuselage 

GRE3 Elastic Rotor Blades 

CSF I Finite Element Structure (Landing Gear and Blade Weight) 
CES! Elastic Stop (Droop and Flap Restraints) 

FRA3 General Aerodynamic Force Module 

STH4 Time History Solver (Integrator) 


The following assumptions apply to the use of DYSCO [Ref. 18]: 
1. All angles are small, which allows approximations. 


2. The relevant physics of a system may be modeled as a set of second-order differ- 
ential equations in the time domain. The equations are of the form: 


MX+CX4+ KX=F 


where X is the vector of the displacements, M, C, and K are coefficient matrices, 
and F is the force vector. i 


3. It is possible to formulate the equations of a system based on the equations of the 
components of the system. 


4. It is possible to compute the state vector (X, X) of each component based on the 
state vector of the system. 


1. Rotor Blade Model 
Since the main concern of this study is the flapping of the rotor blade: ‘he 


flapping equations will be examined in detail. Figure 9 shows the various compornts 


assembled to make the model of the F1-46. 
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Figure 9. DYSCO Model of CH-46 


Figure 10 on page 17 shows the three primary degrees of freedom that a: of 


interest in this study. Ifamiltons’s law of varying action has 
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Figure 10. Degrees of Freedom for Rotor Blade Element. 
Source: DYSCO Theoretical Manual 


been used to derive the equations of motion in a generalized coordinate system which 
allows a direct solution without the consideration of force equilibrium [Ref. 18: p.12]. 
To avoid higher order terms that might complicate the equations of motion, an ordering 
scheme was adopted to determine which terms should be ignored. The following scheme 


was employed in DYSCO: 


Table 1. ORDERING SCHEME USED IN DYSCO. 
ee 












17 


Xp Yu. ZH. Ay Ay a,= O(e) 


where the following are: 


R Rotor blade length, also inertia frame 

u,V,W elastic displacement in the x,y,z directions. 
od Pitch angle of blade element 

Cs principal axes of local cross section of blade 
c chord length of blade element. 

t thickness of blade element 

o local angle of attack of blade element 

e mass centroid offset from elastic axis 

Bic precone angle of the blade 


The references for these variables are shown in Figure 11 on page 19, and Figure 12. 
In Figure IJ, the X,Y,Z coordinates are in inertial frame, R. The R system coordinates 
are the rotor shaft axes when there is no hub motion, while x,, y, s coordinates are fixed 
in the reference frame B, which rotates with respect to R frame at a constant angular 
velocity 2. The coordinates x’, y’,z’ are the coordinate axes of the deflected blade ele- 
ment. 

As mentioned earlier, the Hamilton’s Law of Varying Action is employed to 
derive the equations of motion for the rotor blade. “Hamilton’s Principle”, can be ex- 


pressed as: 
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Figure 11. Blade Coordinate System for DYSCO 
Source: DYSCO Manual, Vol. I 


[ew =i Oran 
tl 


where U is the strain energy, T is the kinetic energy, and 6W 1s the virtual work of the 
external forces. Hamilton’s law is a special case of Hamulton’s Principle, and is [Ref. 
13e9-p: 10) 


| wu- 1) - 5wyar - FT sail =0 


Further developments of strain energy, kinetic energy and the virtual work of 
the external forces are made in Reference 18, and are applied to the blade elastic dis- 
placement through the Raleigh-Ritz method. For this, arbitrary functions for the blade 


displacement can be separated into a sum of products of functions of r and t only: 


Wr) = > OYA) = > WY 
l 

wrt) = > 2)Zlr) = )4Z; 
y 


br) = > POOL) = > OO, 
k 


where Y(r), Z(r), (7) are modal functions and y, z, and @, are generalized coordinates. 
These sums are substituted into the displacement equations for u, v, and w to yield 
equations of motion for all the generalized degrees of freedom.[Ref. 18: p.49] 

These methods are incorporated in the elastic blade module, named CRE3. The 
module allows out-of-plane bending (flapping), in-plane bending (lead-lag), and torsion 
(pitch change)[Ref. 18: p.50]. Mode shape data is also defined, based on the normal 
modes of a non-rotating beam for given boundary conditions[Ref. 19: p.70]. Rotor hub 
degrees of freedom are user defined, and will automatically couple to the fuselage com- 


ponent degrees of freedom through a detailed set of transformation matrices. 
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Figure 12. Hub and Perturbation Rotational Degrees of Freedom. 
Source: DYSCO Manual, Vol. I 
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2. Rotor Blade Aerodynamic Model 

As discovered by Johns in flow visualizations [ Ref: 17] and by Anderson in ve- 
locity measurements [Ref. 6 ], the airflow in the vicinity of the flight deck of a ship is 
extremely turbulent and unsteady. Rhoades [Ref. 5] study is the best effort to obtain 
airflow velocity measurements so far, and his data were used as input to the FRA3 
module. The aerodynamics of the rotor blade are further complicated by the low RPM 
of the rotor blades during rotor engagement, thus introducing Reynolds numbers far 
below those encountered in normal flight. While simple linear aerodynamics (compo- 
nent FRAO) was used for the DYSCO model test runs, a more comprehensive approach 
is needed to accommodate the changing aerodynamics that exist on the ship flight deck. 
This can be accomplished with the FRA3 component of DYSCO. Using this component 
implies that linear aerodynamics does not accurately model the flow over the deck of the 
AOR. Figure 13 on page 23 illustrates the rotor aerodynamic logic for DYSCO. 

All of the steady state aerodynamic coefficients for the rotor blade (C, Cp Cy) 
can be either calculated or looked up from tables external to the DYSCO program. For 
either option, the angle of attack (a) and Mach number (M), based on the two- 
dimensional strip theory, are modified by yaw flow angle to account for the three- 
dimensional aerodynamic effect.[ Ref. 18 : p.52] 

As observed by Rhoades [Ref. 5 ], angles of attack for the upflows from the side 
of the ship and downflows from the forward superstructure are all greater than the low 
Reynolds number stall angles in Reference 20. This indicates that drag forces on the flat 
plate area of the rotor blade are likely to be the largest contributor to forcing the blades 
to such high flapping angles. This view is reinforced by the results of the computer 
simulations and by the blade element program in Appendix A, which shows the lift and 
drag forces as a function of blade position in the rotor disk. The modified drag coefh- 


client 1S 


Con = Co(@rzp; Mey) 


t + 


where agp 1S the angle of attack corrected for compressibility effects. 


B. MATHEMATICAL MODELING OF AIR-WAKE TURBULENCE 
In order to fully understand the modeling methods used for representing the ship 
airwake turbulence, one must understand the use of stochastic control theory. In brief 


terms, the modeling is done by shaping Gaussian white noise through the use of first or 
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Figure 13. Rotor Aerodynamic Logic. 
Source: DYSCO Theoretical Manual 


second order lag functions so that the output closely resembles the spectral densities of 


the velocity components at each point in the wake of the ship. 
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White noise is a theoretical abstraction used to simplify calculations. It is simply a 
random process with an expected value (mean) of zero and an absolutely flat power 


spectrum. 


S(w) = W = constant for all w 


where S(w) is the power spectral density function with respect to frequency w. Because 
the mean square value (p(0)) of any random process is the integral of the spectral density 


over all frequencies, 


(0) =e |” slow) de 


and since white noise has a constant spectral density for all frequencies, white noise 
theoretically has an infinite mean square value. This bothersome feature means that the 
power spectral density does not decrease with increasing frequency. Thus shaping 
functions are needed to represent real-life processes. [Ref. 21] 

The best examples of the modeling of atmospheric turbulence are the Von Karman 
and Dryden models, which are briefly discussed in Reference 1. The models resemble 
free-air turbulence closely enough to be used in the military aircraft design process 
[Ref. 3]. According to the Dryden model the spectrum of the vertical component of 


random wind velocity in turbulent air is 


1+ 3(wT) 


S(w) = oT 
= eT eat y 


where a? is the variance in the z component, and T is the time constant. This spectral 
density function is shown in Figure 14 on page 25. Figure 15 on page 26 illustrates the 
realization of a signal with a Dryden spectrum. 

The two functions being considered here are first and second order lag functions. 


They are 


A 
GAS) = T+ BS and 


K(S + a) 
INS) Sg 
S” + 2¢w,S + w, 
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Figure 14. Dryden spectrum of air turbulence. 


Source: Freidland, B., Control System Design 


where K and A are filter gains, “a” is the location of the minimum phase zero, B and w 
are the frequencies, and ¢ is the system damping coefficient. By adjusting the bandwidth, 
a model representing the measured data can be found. Once the model is formed, the 
transfer functions representing the model may be used to input turbulent disturbances 
into the plant matrices of a flight simulator. To obtain the simplest possible simulator 
model of the turbulent flow, transfer functions must be developed for the u, v, and w 


velocity components. 
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Figure 15. Realization of signal having Dryden spectrum. 
Source: Freidland, B., Control System Design 
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II. CONSTRUCTION OF TUNNEL STRIKE AND AIR-WAKE 
TURBULENCE MODELS 


A. AIR-WAKE TURBULENCE MATH MODEL 

The goal of the air-wake model is to mathematically simulate the turbulence felt by 
a helicopter flying downwind of a ship. The methods used to do this involve finding the 
spectral densities of the air-wake at discrete points along some glidepath, and then de- 
veloping filters to shape random or Gaussian white noise to them. The filter will be of 


the form: 


(5) __K(s—2(1)) (s—2(2))... (s— 2m) 


H(s) = p(s) (s—p(1)) (s— p(2)) ... (s — p(n)) 





Filter design is accomplished here using the Signal Processing Toolbox of MATLAB. 
For this project, an infinite impulse response digital lowpass filter is used, and pre- 
sented in the form given above. The algorithms used by MATLAB take the spectrum 
shape specifications given by the user and produce vectors for the numerator and de- 
nominator, yielding a transfer function in the first canonical form. The software can 
then be programmed to produce random white noise, and pass this noise through the 
filter. The digital domain representation is . 


Y(z) _ 61) + hea en za 
X(z) eae ocean 1)z5 


where a(i) and 6(7) are the transfer function coefficients. 

This filtered white noise is then processed to find the spectral density, using break 
frequencies input by the designer. This effort is directed to providing a transfer function 
for the DD-963 class ship, with data obtained from the Naval Postgraduate School wind 
tunnel measurements being used as the baseline. To obtain the math model, the gains, 
poles and minimum phase zeros are derived for points along helicopter glide paths. 

Appendix D contains the file used to transform the transfer functions to power 
spectra models. The simulated spectra contrast with those offered by Nave [Ref. 10: 
pp.29-32.] and Hanson [Ref. 4: pp.12-29.] in that the curves are smoother, due to 
Healey’s data[Ref. 11] for the spectra containing 66,000 samples, acquired at 2.5 kHz. 
vice 256 samples acquired at 164 Hz for Nave’s FF 1052 model and 131 samples at 164 


Zi 


Hz. for Garnett’s DD-963 model. The lower number of samples and lower sampling rate 
produces a very jagged power spectral density curve, while the higher number of samples 
will give a more accurate variance and mean, and the higher sampling rate will smooth 
out the curve somewhat. Scaling of the bandwidth 1s accomplished by computing the 
ratio between model and prototype size. Garnett used the scale factor of the ship model 
to scale the frequencies of the turbulence up to the prototype. This was justified by the 
use of real-world wind velocities in the wind tunnel. Thus, the 1:50 FF-1052 model's 
frequencies were divided by fifty [Ref. 10: p.23.], and the 1:80 DD-963 model’s fre- 
quencies were divided by eighty [Ref. 9: p.20.]. NPS’s model employed much slower 
tunnel velocities than real-world values, hence the frequencies are scaled using velocities 
and characteristic lengths, in this case the beam of the ship. By comparing the Strouhal 


number for the prototype and wind tunnel model, a scaling factor of fourteen was found. 


B. DYSCO CH-46 TUNNEL STRIKE MODEL 

As mentioned in Chapter II, various modules were used to construct the math model 
of a complete CH-46 helicopter. Rotor blade properties were obtained from the Dy- 
namics Group of Boeing Helicopter Company, and are provided in Appendix B. Ap- 
pendix E contains the information entered into each of the DYSCO modules. To ease 
the complexity of the system, only a two-dimensional rigid body representation was used 
for the fuselage. Parameters such as C.G. location and locations and inclinations of the 
rotor masts were obtained from the CH-46 Maintenance Information Manual [Ref 
2a. 

Another simplification involved eliminating the forward rotor head. Since there 
have been no incidences of tunnel strikes of the forward head, due mainly to its tilt away 
from the fuselage, such a simplification was justified, and lessened the computation time. 

Degrees of freedom deemed important for the aft rotor head are the out-of-plane 
or flapping (OP) degree of freedom, in-plane or lead-lag (IP) degree of freedom and 
torsional or blade pitch (TOR) degree of freedom. All three basic degrees of freedom 
were included to allow for any coupling that may occur between the three. Two modes 
were allowed in the flapping degree of freedom to insure an accurate description of the 
bending modes of the rotor blades. 

To model the flap stops incorporated in the CH-46 to restrict the amount of flapping 
at the hub, the DYSCO module CES1 was used. The module not only imposes re- 
Strictions on the flapping, but also was used to model the lead-lag dampers for each 
blade in the in-plane degree of freedom. Values for the damping, obtained from Boeing 


28 


Helicopter Co., are included in the DYSCO input list in Appendix D. Figure 16 on page 
30 illustrates a sample model using CES]. 

Another constraint was needed for the pitch or torsion degree of freedom. This was 
accomplished by using the simple module CCEO. This module helps model the pitch 
change control rod for each rotor blade. Failure to add this constraint would result in 
a rotor blade that theoretically could twist 360° in torsion. By adding a control rod with 
a very large stiffness value, such twisting is prevented mathematically. 

Linear constraints are added to the fuselage model to simulate landing gear. Using 
CSF1, simple spring damper systems are constructed to model the nose and main land- 
ing gear. Stiffness and damping coefficients are used to ensure that high oscillations do 
not occur that would affect the movement of the rotor blades. These simple spring- 
damper systems are then coupled linearly to the fuselage degree of freedom ZCG, which 
is the CG degree of freedom in the Z direction. A sinusoidal forcing function is added 
to the landing gear system to simulate wave action on the ship, but has been set to zero 
for this study. Future studies could incorporate this forcing function. 

Modeling the air flow from the side of the ship up through the rotor system has been 
difficult. As covered in Chapter II, the module FRA3 will compute aerodynamic forces 
of airflows. A look-up table of low Reynolds’s Number airfoil data was obtained from 
Reference 20, based on a Reynolds number of 660,000. This was used based on the 
computed Reynolds Numbers for each blade at each time step using the algorithm in 
Appendix A, which was based on wind tunnel measurements by Rhoades [Ref. 5], and 
scaled to the prototype. Values of wind-tunnel velocities from Rhoades are given in 
Appendix F. Rhoades’ flow visualization on the windward side also indicated that an- 
gles of the flow ranged from 15° to 25°, relative to the rotor disk. Since these angles are 
well above the stall angle of attack for the airfoil used (12°) the predominant force acting 
on the rotor blade is expected to be the drag force caused by the updraft or downdraft 
of the vertical component of the wind. Appendix G contains the values used in the look 
up table and the inputs for the FRA3 module. 


C. SIMPLE CH-46 BLADE ELEMENT PROGRAM — 

The blade element program in Appendix A is an extremely basic program and was 
developed to describe the flowfield around the rotor disk and over the flight deck of an 
AOR class ship. The program was further modified to compute the aerodynamic forces 
acting on each rotor blade. This program assumes that the blade is a rigid beam, that 


is, no bending is allowed. The algorithm was developed from wind tunnel measurements 
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Figure 16. CESI Model (example) 
Source: DYSCO Users Manual 


made by Rhoades [Ref. 5 ], which are presented in Appendix F. Using the rotational 
acceleration obtained from a videotaped Dynamic Interface wind limitation flight test, 
the varying linear velocities are computed as a function of one-half second time steps for 
three revolutions of the rotor head. These velocities were combined with the flow field 
to obtain the relative wind to the blade element, and the angle of attack at each blade 
element. The following simple calculations illustrate the steps incorporated in the pro- 


gram. 


For ‘Y= 3°, Rotor traveling clockwise, 
Point | at Blade Root 
Point 5 at Mid-Span of Blade 
Point 10 at Blade Tip 
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The curve fit for the rotor RPM acceleration, and subsequent integration to find blade 


position is as follows, 
RPM = 4.2675E—02 + 3.5366¢ — 0.457881? + 2.3410E— 020° 


POS = 0.00712 + 0.029471? — 0.002543r? + 0.0000987* + 0.151325, 


where POS is the number of revolutions completed by that blade. For the 15th time 
step, the time is 7 seconds, which yields a RPM of 10.41, and a blade position of 2.92 
degrees. The horizontal (u) and vertical (w) components of the wind over the flight deck, 
relative to the rotor disk plane, are determined at each position of the rotor blade. For 
this instance, they are found to be 2.0429 m/s and 0.9377 m/s respectively. These values 
are then scaled to represent the 40 knot wind represented by the wind tunnel, which 
yields a horizontal component of 58.6 feet per second and a vertical component of 26.87 
feet per second. 

The circumferential velocities due to rotation are determined by the RPM and the 
position of the blade element from the center of rotation. The following values were 


obtained for time step 15, 

Radius, = weer inches, Vasrerionoi) = 3-O25fps 

Radius, = 148.1 inches, Waprerione’s = 13.45 fps 

Radius,, = 291.64 inches, Voasretionatio = 20.49 fps 

Adjustments are made for the position of the leading edge of rotor blade with respect to 
the direction of the flow field, so that only the components of the flowfield normal to the 


leading edge of each blade element and normal to the chordline of each blade element 
are considered. This is accomplished by 


Uolade = "rotation ~ “wing COS W 


where u,,4, 1s the horizontal wind component normal to the leading edge. For example, 
for the 10th blade element, u,,,,,18 —32.02 fps, indicating that the horizontal component 
is flowing up the trailing edge of the blade element. 
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The inflow angle to the blade element is then determined by computing the 
arctangent of the vertical and horizontal flows at the blade element. Since the CH-46 
rotor blade has a linear geometric twist built in, this twist is also computed and added 
to the inflow angle for each blade element. The program then uses the angle of attack 
to search the look-up table presented in Appendix G and obtain the correct lift and drag 
coefficients each blade element. From this, lift and drag forces are computed with re- 
spect to the relative wind, and finally each are resolved into a vertical component. Re- 
sults of this process are presented in Chapter IV and Appendix A. 


IV. RESULTS 


A. AIR-WAKE TURBULENCE MATH MODEL 

In order to obtain the best physical fit of the computer generated power spectra, 
several orders of transfer functions were attempted. Basing the breakpoints on the slope 
changes of the measured power spectra, a sixth order transfer function was finally used 
to model the measured spectra. The parameters used in the “yulewalk” algorithm in 
“MATLAB?” for the modeling of the DD-963 spectrum at a point twenty-five percent of 
the ship length away from the flight deck, and with the ship pointed thirty degrees to the 


right of the wind are as follows for each component of the velocity in the air-wake: 


m,=[10.0 1.0 —2.0 —32.0 —72.0 —82.0] 
f=[0 3/1000 35/1000 600/1000 999/1000 1] 

m, =[0.0 0.5780 —4.0462 —27.1676 —41.6185 —58.9595] 
f=[0 3/1000 35/1000 600/1000 999/1000 1] 
m,=[0.0 1.0 —-3.0 —27.0 —45.0 —65.0] 


f =[0 3/1000 35/1000 600/1000 999/1000 1] 


where 7m iS a vector containing the desired magnitude response at the points specified in 
f(breakpoints) for each of the components (u, v, and w) and f is a vector of frequency 
points, specified in the range between 0 and 1, where 1.0 corresponds to half the sample 
frequency (the Nyquist frequency). The transfer function coefficients (see p. 28) from 


those models are: 


b, = [22.6588 —13.0029 —9.0971 —4.4572 —0.1487 6.3585 —0.4214] 
a@,=[1.0 0.5978 —0.1317 —0.3804 —0.3665 —0.0400 0.0205] 
b, =[19.0600 —6.0289 —9.0813 —2.9458 —1.5902 2.8828 0.4594] 


a,=[1.0 0.5388 —0.1256 —0.2584 —0.2770 —0.0516 0.0040] 


o5 


b,, = [18.5670 —3.7482 —5.4913 —8.4648 —4.1090 5.0879 0.8924] 
a, =[1.0 0.7695 0.2652 —0.2462 —0.4730 —0.1138 0.0031] . 


These are placed in the following canonical form of a transfer function, with Y(z) as the 
output and U(z) as the input. 


Y(z) (1) + O(2)z7) +... + Bt Dz” 


U(z) 1 + a2\z 09 + ey 





Having modeled the power spectrum of the wind tunnel measurements, the model 
was then scaled to the prototype using a Strouhal scaling factor of fourteen. This was 
further reduced to a pole, zero, gain format transfer function, using conversion algo- 
rithms in MATLAB. This format is useful in determining the inputs needed for digital 
simulation. The measurements used for this simulation were the u, v, and w components 
of the velocities of the air-wake at a point twenty-five percent of the ship length away 
from the touch-down point on the flight deck. These measurements are shown in 
Figure 17 on page 35. The transfer coefficients, when used to filter random generated 
noise, produce the spectra shown in the following figures. The solid line in each spectra 
is the MATLAB generated model, while the dashed line power spectra is the measured 
spectrum by Healey [Ref. 11 ]. Both the wind tunnel and Strouhal scaled models are 
presented. 

When comparing the measured and simulated spectra, the smoother nature of the 
former is noted. This could be attributed to the very large size of the samples of meas- 
ured data, which were then averaged in groups of sixty-four, while the simulated spec- 
trum curve was formed with a random signal generator that had no way to control the 
variance of the signals. The larger variance in the simulated spectra could be a source 
of higher than normal turbulence, but should be evaluated by a real-time simulator flight 
prior to adjustments. The MATLAB predicted spectra are much smoother than those 
obtained by Fortenbaugh and Hanson. For example, the power spectral density curves 
derived by Hanson for yaw angle of 330° shows very large variances. The very wide 
variance presented by Hanson could be due to his small data size and small sample rate 
[Ref. 4: p.30.]. Figure 24 on page 42 is presented only to illustrate the greater variance 
(more “spikey”) at the lower frequencies in the Hanson and Fortenbaugh models. In the 
NPS wind tunnel, frequencies lower than three radians per second were not measured, 


but the overall smoothness Healey’s measured spectrum indicates a lower variance. The 
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Figure 17. Wind Tunnel Spectra - U, V, and W Components 
algorithin used by Ianson allowed for changing the variance of the air wake component 


velocities, thus permittig him to reduce the variance by the sixty to seventy percent 
needed to comply with pilot’s reports of the simulator trials [Ref. 4: p.30.]. 
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Figure 18. Filtered Noise for U Component Simulation. 
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S(f), VY Component 


V-Component Spectral Density, DD-963, 30 Degrees Yaw 
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Figure 19. Filtered Noise for V Component Simulation. 
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S(f), W Component 


W-Component Spectral Density, DD-963, 30 Degrees Yaw 
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Figure 20. Filtered Noise for W Component Simulation. 
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S(f), U Component, Prototype 


U-Component Spectral Density, DD-963, 30 Degrees Yaw 
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Figure 21. Filtered Noise for U Component Prototype. 
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S(f), V Component, Prototype 


V-Component Spectral Density, DD-963, 30 Degrees Yaw 
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Figure 22. Filtered Noise for VW Component Prototype. 
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W-Component Spectral Density, DD-963, 30 Degrees Yaw 
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Figure 23. Filtered Noise for W Component Prototype. 
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Figure 24. Comparison of Hanson’s and Fortenbaugh’s Models. 


The high variance values found in the MATLAB simulation could also be reduced 
by increasing the number of points used in the random signal generator. For this model, 
4000 points were used. Larger sizes of white notse sets could not be generated due to 
memory linutations with the PC version of MATLAB. Higher functioning versions of 
MATLAB, on systems with much larger memory could produce sets one or two orders 
of magnitude higher, thus reducing the problem of higher variance in the simulated 
spectra. To test Chis assumption, the same algonthm was run en PRO-MATLAB, in- 
stalled on the VAN Muini-2000 using VAN VMS. Since the memory available wes con- 
siderably greater than that ofa PC, the number of white noise points was increased by 
an order of magnitude, ta 25,000. The power spectra obtained were nwch smoother, and 
gave an mproved representation of the wind tunnel flowfteld. This approach also con- 


firms Tlanson’s theory that the frequency used to generate the white noise signel will 


C 
ae 


attect the quality of the simulation (Ref. 4: p.11j]. By varnving the number of discrete 


White noise points generated, i may be possible to elinunate the need to interpolate the 


white noise prior to filtering when using Hanson’s model. The improved computer 
simulation for the full scale prototype is shown in Figure 25. 

The poles, zeros and gains obtained from MATLAB for the wind tunnel model are 
given in Table 2 through Table 4. 

Since the random signal changes with each computer simulation, the portrayed 
spectral densities can be assumed to change somewhat with each real-time run. Further 
study in conjunction with real-time simulator flights are needed to further refine this 
approach to air-wake simulation. Results of filters designed for the complete glidepath 
for a 330° yaw angle are presented in Appendix H. 


Table 2. POLES, ZEROS AND GAINS FOR U COMPONENT OF THE WIND 
TUNNEL MODEL 
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Table 3. POLES, ZEROS AND GAINS FOR V COMPONENT OF THE WIND 
TUNNEL MODEL 
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B. DYSCO CH-46 TUNNEL STRIKE MODEL 
After several attempts to accurately duplicate the velocity field around the whole 
rotor disk, it was determined that a separate flowfield would be considered for each rotor 
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Comparison of Two Models with Different Variances 


10-2 


10° 


S(f) 


10+ 





10° 
10! 10° 10! 102 


Frequency (Hz) 


Figure 25. Wind Tunnel U-Spectra Simulation Using PRO-MATLAB. 


blade. This assumes that the behavior of the rotor blades are independent of each other, 


without any interaction via the hub, since the rotor head ts fully articulated. It also as- 


Table 4. POLES, ZEROS AND GAINS FOR W COMPONENT OF THE WIND 
TUNNEL MODEL 
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sumes that the airflow velocity conditions at a particular point in time is constant for the 
length of the rotor blade. Since the conditions of the velocity field were dependent on 
the rotor blade position, a history of rotor blade RPM vs. time was required. It was 
obtained from videotape of a Dynamic Interface flight test conducted aboard an AOR 
class ship. The time history is illustrated in Figure 26 on page 46. 
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Figure 26. Rotor RPM Time History for Rotor Engagement 


By combining this time history with the flowfield measured by Rhoades [Ref. 5], the 
velocity components relative to the rotor disk for each time step were obtained. This 
fixes the field in space, and uses only one Reynolds number for the range of velocities 
each rotor blade will “see.” Actual relative velocities (combined effects of mean airflow 
and rotor blade velocity) are computed by DYSCO. The algorithm for mean velocity 
components and blade position is based on the velocity measurements obtained for the 
AOR wind-tunnel model [Ref. 5: p.71.]. The entire program used to determine the ve- 
locities, and scale them to the 40 knots simulated in the wind tunnel is listed in Appendix 
A, along with the results for each rotor blade. 

These velocity components were then input to the DYSCO CH-46 model by 
changing the velocity inputs in the FRA3 module for each time step. The wind values 


were input using downward flows as positive (+) and upflows as negative (-), consistent 
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with the DYSCO users manual [{Ref. 18: p.56.]. When combined with lift, drag and 
moment coefficients for a NACA 0012 airfoil, the Runge-Kutta integration in DYSCO 
could not converge to a solution. After several attempts to correct the divergence, use 


of DYSCO for the tunnel strike model was abandoned. 


C. SIMPLE CH-46 BLADE ELEMENT PROGRAM 
Using the methods outlined in Chapter IV, the total vertical forces for each rotor 
blade where computed with respect to the rotor blade position. These forces are pre- 


sented in Figure 27, Figure 28, and I'rgure 29. 
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Figure 27. Total Vertical Forces on Blade One. 


These results indicate that for blade one, there are large upward forces at the 180° 
postion and considerably smaller downward forces (less than 50 Ibs.) at the 360° post- 
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tion. Although this program does not include centrifugal forces on the rotor blade, such 
forces would limit the affect of the large upward forces on blade displacement. Blade 
two exhibits similiar behavior, with large (150 Ibs.) upward forces at the the 180 ° posi- 
tion and low downward forces at the 360° position. Blade three exhibits forces opposite 
that of one and two, showing a large upward force near the 360° position, and a small 
downward force near the 220° position. The distribution of forces agree with the be- 
havior of the blade displacements observed in tunnel strike videos. The shift in forces 
on blade three are due to the increased rotational velocities when the strong upflow is 
encountered. Beyond the second revolution, the linear velocities due to the rotation of 
the blades prevent the net forces from acting in the downward direction. Although the 
aerodynamic forces appear to be increasing past two revolutions, the resultant force on 
each blade would move away from the vertical as the centrifugal forces on the rotor 
blade increased with RPM. The aerodynamic forces correspond with known behavior 
of the rotor blades in low RPM environments, where the blades tend to rise or “sail” at 
the 180° position, and dip down when the blade is over the fuselage of the aircraft 
(360°). 
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Figure 28. Total Vertical Forces on Blade Two. 
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Vertical Blade Three Forces 
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Figure 29. Total Vertical Forces on Blade Three. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. AIR-WAKE TURBULENCE MATH MODEL 

Initial results show good agreement between predictions of the sixth order math 
model and the measured data from the wind tunnel model. Higher order transfer func- 
tions appear to model the measured turbulence with greater accuracy. Such a high order 
transfer function would not present computational difficulties in a real-time simulation, 
since the power spectra could be calculated prior to the running of the helicopter simu- 
lator, and the results stored for use as the helicopter passed through the points simulated 
by the math model. Since each point would have only thirty-nine variables (eighteen 
poles, eighteen zeros and three gains) for the components of the turbulence, the prior 
computation would not be time consuming nor utilize a large area for data storage. 

This approach requires the measurement of the airwake at numerous points along 
the prescribed approach paths for each type of helicopter the DD-963 class of ship. 
Measurements would also have to be made with the wind at various angles to the ship 
heading. For those measurements already completed, filters must be designed to model 
the power spectra for each component. As noted in Chapter IV, comparison of the re- 
sults with Hanson’s first order model showed a much smoother curve (with less vari- 
ance). This results in an overall lower variance than Hanson’s model, eliminating the 
need to reduce the variance of the turbulence velocity components to simulate real-time 
conditions. The scaling of the frequencies to match the frequencies of the prototype will 
not affect the coefficients of the transfer functions. The output is simply scaled to match 
the lower frequencies of the prototype. Thus, the poles and zeros of the transfer func- 
tions will not change when applied to the prototype. 

The smoothness of the computer model is enhanced with a smaller variance, ac- 
complished by using a larger set of random signals to simulate the white noise. As the 
number of white noise points increases, the variance approaches zero, thus reducing the 
scatter in the PRO-MATLAB simulation that occurred in the PC simulation. Needing 


further study are: 
Computational compatibility with existing real-time simulators. 
Measurement and modeling of additional approach paths to the ship’s flight deck. 


Pilot evaluation of the turbulence model on a real-time simulator. 


ai 


Since actual frequencies over 1.0 Hz. are not considered important for most manned 
air-vehicles, the response of the filters at those higher frequencies are considered to be 


negligible. 


B. CH-46 TUNNEL STRIKE MODELS 

In this preliminary analysis, an attempt was made to “produce” a blade strike by 
exposing the helicopter blade to the mean flow, as measured by a hot-wire anemometer 
over the flight deck of an AOR model. Due to the extreme divergence occurring during 
the Runge-Kutta integration scheme in DYSCO, a solution was not obtained. Further 
study is needed to investigate if the DYSCO routines are suitable for such low Reynolds 
Number applications. The results of the wind tunnel measurements should also be used 
in more sophisticated blade analysis programs, such as RACAP and CAMRAD. For 
the Blade Element Program, progress was made in identifying those areas of upward and 
downward aerodynamic forces on the rotor blades. The program has shown that meas- 
urements of flowfields over model ships in wind tunnels can be applied to estimating the 
distribution of mean aerodynamic forces on each rotor blade with respect to blade po- 
sition. Further development of this program should include its validation, incorporation 
of rotating beam theory, accounting for the blade mass and bending characteristics along 
with rotor blade flapping. Further measurements of flowfields over various classes of 
ships should also be made at various yaw angles to analytically investigate those areas 


of high aerodynamic forces present during rotor engagement and disengagement. 
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APPENDIX A. CH-46 BLADE ELEMENT PROGRAM 
PROGRAM ELEMENT1 


Tete dedevededevevede se seve see TRI RRR RRR RK RI RR RSS Se 


* 
* 
* 
* 
¥ 


PROGRAM TO DETERMINE VELOCITY FIELD * 
OF AOR FLIGHT DECK AT 90 DEGREES YAW * 
FOR GIVEN ROTOR ENGAGEMENT ACCELL * 
AND FIND LIFT AND DRAG FORCES * 


ON THE ROTOR BLADE USING BLADE ELEMENTS* 


tetetetetetetevedetetetevedete tetete tek fete Lek FF KLEE RIECK ERERERE 


DIMENSION TIME( 41) ,U(41) ,W(41) ,POS(41) ,REV(41) , PHI(41),RPM(41), 
+VROT( 41,10) ,UM(41) ,WM(41) ,AOA( 41,10) ,AOAL( 101) ,CLL(101) ,CL( 41,10), 
+CDL(101),CD(41,10) ,AOAT( 41,10) ,DRAG(41,10) ,RSEC(10) ,TWIST(10), 
+DBLADE( 41) , AOAAVG( 41) ,AOAD( 101) , DRAGV(41,10) , DBLADV( 41) , FORCE(41) 

REAL MAG1(41,10),LIFT(41,10) , LBLADE(41) , LIFTV(41, 10), LBLADV(41) 

DOUBLE PRECISION MU, REN1(41,10),RENAVG(41) 

OPEN( UNIT=9 , FILE='FLOW1. DAT’ , STATUS='old' ) 

OPEN( UNIT=10, FILE='TIMESTEPS. DAT' , STATUS=' OLD' ) 

OPEN( UNIT=15 , FILE=' LCOEFFE. DAT' , STATUS='OLD' ) 

OPEN( UNIT=17 ,FILE=' DCOEFFE. DAT’ , STATUS=" OLD ) 

OPEN( UNIT=16 , FILE='SECTIONS. DAT , STATUS=' OLD' ) 

OPEN( UNIT=18 , FILE='GRAPH1. DAT' , STATUS=' OLD' ) 


REV=NO. OF REVOLUTIONS 

POS=ANGLE OF REFERENCE BLADE IN RADIANS 
=VELOCITY COMPONENT TANGENTIAL TO ROTOR DISK 
=VELOCITY COMPONENT PERPINDICULAR TO ROTOR DISK 


READ(10,*)(TIME(I), I=1,41) 
READ(15,*)(AOAL( J) ,CLL(J), J=1,101) 
READ(17,*)CAOAD( J) ,CDL(J), J=1,101) 


COMPUTE BLADE POSITION 


DO 100 I=1,41 

REV(I)=. 0071*TIME(I)+0. 029470*( TIME(I) )**2. 0-0. 002543*(TIME(I)) 
+**3. 0+0. 000098*( TIME(T) )**4. 0+0. 151325 

RPM(1)=4. 2675E-02+3. 5366*TIME(1)-0. 45788*TIME(I)**2. 0+2. 3455E-02 
+*TIME(I)**3. 0 


CHANGE REVOLUTIONS TO RADIANS 


POS(I)=REV(1)*2. 0*3. 141593 
PHI(1I)=REV(1)*360. 00 


COMPUTE THE TWIST ANGLE AND RADIUS OF EACH BLADE SECTION 


DO 25 L=1,10 
RSEC(L)=(L*28. 71-28. 71/2. 0)+18.9 
TWIST(L)=(4. 353-0. 02779*RSEC(L)) 

CONTINUE 
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C 
C COMPUTE THE ROTATIONAL VELOCITY AT EACH BLADE ELEMENT 
C 
DO 27 J=1,10 
DO 28 N=1,41 
VROT(N,J)=RPM(N)*2. 0*3. 1416/60*(RSEC(J)/12. 0) 
28 CONTINUE 
27 CONTINUE 
C 
C COMPUTE VELOCITIES HORIZONTAL AND PERPINDICULAR TO ROTOR DISK 
C 
20 if((PHI(1I). GT. 0.0). and. (PHI(I). le. 90.0))goto 30 
if((phi(i). gt. 90. 0). and. (phi(i). le. 180.0))goto 40 
if((phi(1). et. 180. 0), anda(phica > lem27 ONO) coco 50 
if((phi(i). gt. 270.0). and. CPHI( jo LE 56050) }GaTomco 
if(phi(i). gt. 360. 0)phi(i)=phi(i)-360. 0 
goto 20 


50) UC I)=2. 05+(phiCi)/90. 0)*(1. 83-2. 05) 
WC I)=0. 944+( phi(Ci)/90. 0)*(0. 87-0. 94) 
goto 70 


40 UCI )=1. 83+( (phd (i) -90m0) /90.0)*( 1560-8 §3) 
WC I)=0. 87+( (phi(i)-90. 0)/90. 0)*(0. 84-0. 87) 
goto 70 


50 UC 1)=1. 80+( (phi(i)-180. 0) 790. 0)*( 1261-1. 80) 
WC I)=0. 84+((phi(i)-180. 0)/90. 0)*(0. 53-0. 84) 
goto 70 


60 UC I)=1. 61+( (phiCi)-270. 0)/90. 0)*(€2. 05-1. 61) 
WC1I)=0. 53+(( phiCi)-270. 0)/90. 0)*( 0. 94-0. 53) 
70 continue 
100 continue 
va=(2. 327**2+., 37**2)**O, 5 
do 200 k=1,41 
C CONVERT WINDTUNNEL COMPONENTS TO SCALED IN FPS 
u( K)=(UCK)/va)*(40*1. 688) 
WCK)=(WC(K) /va)*(40*1. 688) 
200 continue 
C 
C COMPUTE THE REYNOLDS NUMBER FOR EACH TIME STEP 
C 
DO 610 N=1,10 
DO 300 I[=1,41 
TERM1=VROTCI ,N)-UCI)*COS(POS(I)) 
MAG1C I ,N)=(TERM1**2.0 + WCI)**2. 0)**0. 5 
REN1( I ,N)=6400* 1. 6*MAG1(I ,N) 


300 CONTINUE 
610 CONTINUE 
C 


C DETERMINE THE ANGLE OF ATTACK FOR EACH TIME STEP 
C 

DO 620 N=1,10 

DO 600 I=1,41 

TERM1=VROT(I ,N)-UCI)*COS(POS(T)) 
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AOAT( I,N)=(ATAN(W( I) /(TERM1)) )*( 180. 0/3. 1416)+TWIST(I) 
IF(AOAT( I,N). LT. 0)AOAT(1I,N)=AOAT(I,N)+180. 0 
DO 700 J=1,101 
IF(AOAT(I,N). EQ. AOAL(J))GOTO 900 
IF(AOAT(I,N). LT. (AOAL(J)+0. 5). AND. (AOAT(I,N). GT. AOAL(J))) 
+GOTO 910 
IF(AOAT(I,N). GE. (AOAL(J)+0. 5). AND. (AOAT(I,N). LT. AOAL( J+1)) ) 
+GOTO 920 
GOTO 700 
900 AOA( I ,N)=AOAT(I,N) 
CL(1I,N)=CLL(J) 
GOTO 2000 
910 AOA(I ,N)=AOAL( J) 
Clie N)=CLL( I) 
GOTO 2000 
920 AOA( I ,N)=AOAL( J+1) 
CL(1I,N)=CLL(J+1) 
GOTO 2000 
700 CONTINUE 
2000 DO 750 J=1,101 
IF(AOAT(I,N). EQ. AOAD(J))GOTO 901 
IF(AOAT(I,N). LT. (AOAD(J)+0. 5). AND. (AOAT(I,N). GT. AOAD(J))) 


+GOTO 911 
IF(AOAT(1I,N). GE. (AOAD(J)+0. 5). AND. (AOAT(I,N). LT. AOAD( J+1))) 
+GOTO 921 
GOTO 750 
C 901 AOA(I,N)=AOAT(I,N) 
901 CD(1I,N)=CDL(J) 
GOTO 600 
CR AOA(I,N)=AOAD(J) 
911 CD(1I,N)=CDL(J) 
GOTO 600 
C 921 AOA( I ,N)=AOAD(J+1) 
921 CD(I ,N)=CDL( J+1) 
GOTO 600 
750 CONTINUE 
600 CONTINUE 
620 CONTINUE 
C 
C DETERMINE LIFT AND DRAG FORCES 
C 


DO 660 J=1,10 
DO 650 N=1,41 
TERM1=VROT(N, J) -U(N)*COS(POS(N)) 
MAG1(N,J)=(TERM1**2.0 + W(I)**2. 0)**0. 5 
LIFT(N,J)=CL(N,J)*0. 5*0. 002378*(MAG1(N,J)**2. 0)*(28. 71/12. 0)*1. 6 
DRAG(N,J)=CD(N,J)*0. 5*0. 002378%(MAGI1(N,J)**2. 0)*(28. 71/12. 0)*1. 6 
LIFTV(N,J)=LIFT(N,J)*ABS(COS( AOA(N,J)*3. 1416/180) ) 
DRAGV(N,J)=DRAG(N,J)*SIN(AOA(N,J)*3. 1416/180) 

650 CONTINUE 

660 CONTINUE 

G 

C COMPUTE SUM OF LIFT AND DRAG FOR BLADE ELEMENTS, AVG AOA AND REN 

c LBLADE(1)=0. 0 
DBLADE(1)=0. 0 
AOAAVG(1)=0. 0 


> 


RENAVG(1)=0. 0 
LBLADV(1)=0. 0 
DBLADV( 1)=0. 0 
DO 1500 N=1,41 
DO 1000 I=1,10 
LBLADE(N)=LBLADE(N)+LIFT(N, 1) 
DBLADE(N)=DBLADE(N)+DRAG(N, I) 
LBLADV(N)=LBLADV(N)+LIFTV(N, I) 
DBLADV(N)=DBLADV(N)+DRAGV(N, I) 
AOAAVG(N)=(AOAAVG(N)+A0A(N,1)) 
RENAVG(N)=(RENAVG(N)+REN1(N,1))/I 
1000 CONTINUE 
AOAAVG(N)=AOAAVG(N) /10 
FORCE(N)=LBLADV(N)+DBLADV(N) 
1500 CONTINUE 
WRITE(9,10) 
10 FORMAT(3X, 'TIME(SEC)' ,2X, 'POSITION' ,2X,'RPM',3X,'AOA AVG' ,6X, 
+'REN AVG' ,6X,'LIFT' ,3X, 'DRAG') 
WRITE(9,15)(TIME(I),PHI(1),RPM(I) ,AOAAVG(1) ,RENAVG(I) , LBLADV(I), 
+DBLADV(I), I=1,41) 
15 FORMAT(F8. 2,2X,F8. 2,2X,F7.2,3X,F7. 2,3X,D10. 4, 4X, 2F8. 2) 


WRITE(16,*)(' BLADE POSITION ANGLE OF ATTACK’ ) 
WRITE(16,1600)(PHI(N),(AOA(N,I), I=1,10), N=1,41) 
WRITE(16,*)('BLADE POSITION REYNOLDS NUMBER’ ) 
WRITE( 16,1700)(PHICN),(REN1(N,I), I=1,10), N=1,41) 
WRITE(16,*)('BLADE POSITION RELATIVE VEL. IN FPS. ") 
WRITE( 16,1800)(PHIC(N),(MAG1(N,I), I=1,10), N=1,41) 
WRITE(16,*)('BLADE POSITION SECTION LIFT (LBS. )') 
WRITE(16,1800)(PHI(N),(LIFT(N,I), I=1,10), N=1,41) 
WRITE(16,*)('BLADE POSITION SECTION DRAG (LBS. )') 
WRITE( 16,1800)(PHI(N) ,(DRAG(N,I), I=1,10), N=1,41) 
WRITE(16,*)('BLADE POSITION SECTION VERTICAL LIFT (LBS. )') 
WRITE(16,1800)(PHI(N),(LIFTV(N,I), I=1,10), N=1,41) 
WRITE(16,*)('BLADE POSITION SECTION VERTICAL DRAG (LBS. )') 
WRITE(16,1800)(PHI(N),(DRAGV(N,I), I=1,10), N=1,41) 
WRITE(16,*)(' BLADE POSITION LIFT COEFFICIENT' ) 
WRITE( 16,1800)(PHI(N),(CL(N,I), I=1,10), N=1,41) 

WRITE(16,*)( 'BLADE POSITION DRAG COEFFICIENT’ ) 


WRITE(16,1800)(PHI(N),(CD(N,I), I=1,10), N=1,41) 
WRITE(18,1900)(POS(N),FORCE(N), N=1,41) 

1600 FORMAT(F6. 2,2X,10F6. 1) 

1700 FORMAT(F6. 2,2X,10D9. 3) 

1800 FORMAT(F6. 2,2X10F7. 1) 

1900 FORMAT(2X,F6. 2,2X,F9. 3) 
END 
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AERODYNAMIC VALUES OF BLADE ONE 


TIME POSITION RPM AOA AVG REN AVG LIFT DRAG 
(SEC) (DEG) (DEG) (LBS) (LBS) 
. 00 54.48 04 146. 16 -4723D+05 SOO oe 20.00 
0 56-30 1.70 139352 -4091D+05 =20.5) Zhe 
200 66:76 3.14 L262 woo ecu TOS eo.12 eo 2 
oO Doe27 4.40 100. 80 - 2209 N+t05 woul, 290 
. 00 Dee 2.) 5.47 65.16 . 3584D+05 3.96 12. 34 
0) 114. 25 O59 39.96 - 51830705 39.40 28. 70 
. 00 Vopr? 7 ries, Laie 09790705 $0.96 Pie 
20 159. 43 1202 21.24 . 8348D+05 128633 foa2t 
. 00 184. 89 8. 36 19.44 ; 62590705 145. 38 16.67, 
0 ZA1-.86 8.62 TSac0 - 7 OED EOS P26 ir 6 
. 00 240. 10 Gaz 2 ie , 28200705 Dio. 00) 107 
70 269. 43 Ino5 a2a20 - SoZ0DT05 9. 69 Bi? 6 
700 Zoe 2 9.84 118.44 . 2204D+05 =5. 61 S213 
0) So) 015 (he, LOS 145. 44 PS UZ 0> -42.47 i Ware) 2 
. 00 Zeoe 10.41 149. 40 »4853D+05 zoo. 70 28. 08 
DU Sos 10271 140.04 157 Ds O5 “302.55 24. 80 
. 00 C209 11.04 98. 28 O10 D0 5 =e G7 6.08 
0 104. 66 ie4 3 43.56 ~ 0200705 37.66 age) G 
. 00 140.91 11. 88 Zone - 8580D+05 122. 64 Slee Xe! 
. 20 178. 69 E2243 21560 S995 m Ua S 175.47 26.68 
. 00 ae 2G 13206 2a - 8584D+05 122.47 Oe 
0 ZOO 3 135202 95.25 soLge2ut0> 99,17 ee 0 
. 00 304. 34 14. 76 110. 16 - 2472D+05 =5.-01 7.20 
. 90 oo0..66 15. 83 143. 64 72007 Da.0> -46. 29 22.96 
. 00 42.49 1708 122.04 . 3014D+05 =>.) > 14. 38 
> 0 97.40 Poe 2 41.76 . 6787D+05 DO. 29.04 
200 156, 95 20517 ZO 26 21 B70D506 Zoo. 03 CoO) 
0 Zod 9 1 Pe) 19.44 . 1086D+06 ZUG. 9 1 17. 46 
. 00 292590 24.17 03502 . 9235D+05 1932 11. 34 
0 Oe 2 26. 56 120396 2ol2GDs05 S23 16. 43 
. 00 96. 19 29.23 33. 84 : 95 /2070> 95,09 Oued 
20 IOS 32.20 15. 84 . 1536D+06 465. 23 17.42 
. OO 293.03 35. 48 47,46 . 8278D+05 64. 69 11. 68 
0 47.51 ee O 68. 76 1/6920105 25.00 ZoOy 
oO 172. 86 43.07 14. 40 . 1862D+06 674. 24 14. 85 
20 SihO7 D 47.41 D200 . 1013D+06 Oe 5 14. 42 
. 00 102.34 52.14 20. 88 . 1666D+06 “107 12 20762 
50 268. 81 dey 15. 48 . 1682D+06 343.74 4.38 
. 00 91.41 G27a2 Z 160 . 1860D+06 503. 34 Toe 12 
50 29143 68. 81 23.04 . 1801D+06 O70 G.o2 oe 7 1 
700 150. 24 Ta 20 Toe 2 06 1030. 74 | 


Oe ee ee ee ee ee ee ee ee 
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AERODYNAMIC VALUES FOR BLADE TWO 


TIME 
(SE 


C) 


. 00 


bt et ft pt pt tt pt pt ps ps ps 
OO WOAONNANUNEEWWNHNHPrPHrOOWUOMDANNADUNFHEWWNHNHHE 


NO 
© 


a0 
00 
50 
00 
50 
00 


220 
. 00 
530 
. 00 
30 
. 00 
00 
. 00 
. 50 
. 00 
250 
. 00 
70 
200 
20 
. 00 
2510 
. 00 
3010 
. 00 
20 
. 00 
50 
. 00 
5S). 
. 00 
50 
. 00 
wo0 
. 00 
0 
200 
=o) 
. 00 


POSITION 


(DEG) 


174. 
Sie 
186. 
199: 
215). 
234, 
25>). 
249. 
304. 
Bole 
2210 
29. 
hese 
90: 
Zi: 
SD: 
hoo. 
224. 
260. 
258. 
338. 
ZU. 
64, 
tel 
162. 
207. 
ace 
341. 
a2. 
130. 
216. 
o1G. 
Dae 
167. 
22. 
70. 
HiGapse 
20: 
ZA 
S41. 
270. 


48 
29 
76 
27 
a7 
25 
aa 
43 
89 
86 


43 
72 
89 
92 
82 
69 
66 
oy 
69 
28 
03 
34 
66 
49 
39 
98 
a1 
90 
a2 
19 
18 
63 
a 
85 
WS. 
34 
81 
41 
43 
24 


RPM AOA AVG 


(DEG) 

04 28. 80 
ie 25.60 
3.14 25202 
4.40 24. 48 
5.47 24. 48 
639 Zee 
Foe 37.600 
7. 82 78. 84 
8. 36 127. 80 
8. 82 144. 00 
Gree 1S 0merZ 
Wea) 144. 00 
9.84 118. 08 
10s 61.20 
10.41 31. 68 
TOs 23.76 
11.04 21. 60 
11.43 23,76 
11.85 S525 
12.43 108.72 
13. 08 144. 36 
137765 143. 28 
14. 76 9750 
iSa05 34. 20 
17208 20. 88 
Tere 20. 16 
20.17 45.36 
22.05 128. 16 
24.17 90. 36 
20. 56 22582 
29.23 L/22e 
Jee) 72. 36 
35.48 68.04 
39-10 15. 48 
43.07 38. 88 
47.41 39°96 
52.14 12. 60 
Sie 61720 
62. 82 Divo 
68. 81 39.24 
1520 i2oU 
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REN AVG 


, 6566D+05 
. 7021D+05 
. 7297D+05 
. 7237D+05 
. 6685D+05 
.5549D+05 
, 3913D+05 
, 2340D+05 
, 2468D+05 
, 4042D+05 
.5154D+05 
.4219D+05 
, 2975D+05 
.4180D+05 
. 6901D+05 
. 9010D+05 
. 9392D+05 
. 7672D+05 
. 4682D+05 
, 2340D+05 
. 3537D+05 
. 3771D+05 
, 3340D+05 
. 7391D+05 
. 1101D+06 
. 1017D+06 
.5512D+05 
, 2857D+05 
. 4357D+05 
, 1196D+06 
. 1329D+06 
.5971D+05 
.7175D+05 
. 1739D+06 
. 1049D+06 
, 1207D+06 
, 1944D+06 
, 1129D+06 
, 2324D+06 
. 1627D+06 
, 2188D+06 


LIFT 
(LBS) 


104. 
114. 
114. 
105. 
84. 
54. 
Z0n 
og 
=) 1 
-49, 
iol 
-49, 
=r 
w2o 
7O- 
140. 
158. 
96: 
26. 
=o. 
“42. 
-45. 
aoe 
ae 
20 
Pe 
Zhe 
22's 
43 
ZAG. 
D206. 
1. 
oe 
580. 
126. 
147. 
601. 
So. 
709. 
304. 
468. 


29 
81 
27 
65 
55 
05 
45 


Zu 
ie 
35 
70 
af 


00 
78 
92 
20 
1 
64 
07 
74 
43 
55 
37 
84 
25 
98 


14 
68 
78 
48 
68 
35 
80 
35 
31 
20 
a2 
15 


DRAG 
(LBS) 


3/535 
41.10 
29. 86 
235164 
19200 
16.54 
12305 


14.59 
28. 66 
29228 
23.922 
12558 
15222 
31. 56 
29108 
24231. 
Toe 
13565 


21-26 
24. 44 


32.59 
205097, 
18.77 
10255 
15.32 
15. 2 
28.44 
16367 
14. 64 
yA S10) 
Toacu 
11.24 
IES) 


24. 21 


18. 68 


C. AERODYNAMIC VALUES FOR BLADE THREE 


TIME POSITION RPM AOA AVG REN AVG LIFT DRAG 
(SEC) (DEG) (DEG) (LBS) (LBS) 
200 294.47 04 144. 00 .3173D+05 =lO07 2 10. 16 
0 298. 29 70 142. 20 3 1SZDT05 -18.04 12. 36 
1. 00 306. 76 3.14 145.44 - 3500D+05 -26. 96 14.72 
DO aloe? 4.40 149. 76 -4189D+05 =42. 56 T7226 
2.00 Soo 20 5.47 Lod Z . 9021D+05 =65..19 22. 04 
2750 354. 25 6.39 1p aoe - 5674D+05 =B2,05.. 28562 
3. 00 Le? 7 75.16 149. 04 . 5304D+05 =] 52109 32.24 
5250 39.43 Vos 139. 68 . 4059D+05 -40. 20 30.54 
4.00 64. 89 8. 36 111. 60 LoD ODGOD = eco 103 
4.50 91786 B. Oe 59. 76 . 3984D+05 7.94 Lol 
25 Oe) 120. 10 cal 34. 20 . 6363D+05 64.55 32.42 
20 149. 43 oo 24. 48 . 8409D+05 123209 27.64 
6. 00 oe 2 9. 84 21. 60 92590505 154. 35 24.37 
6.50 210269 1G. 13 212916 . 8228D+05 118. 10 19213 
7.00 242. 92 10.41 28. 08 . 5902D+05 Dene ios. 
0 eee lO)e 2a 61.56 -3151D+05 6.54 Volo 
8. 00 309. 69 11. 04 129. 60 . 2432D+05 -13.14 TZ16 
8.50 344. 66 11. 43 148. 32 -4149D+05 -54. 32 22208 
9.00 PAOK SI 11. 88 145.44 - 4099D+05 7 > l30 22709 
g. '0 58. 68 12. 43 113. 40 . 2980D+05 =aeoe 10.41 
10. 00 96.27 13. 08 47.52 - 5475D+05 505 29 26,12 
1025 0 140. 02 13.65 2520 . 9065D+05 Weis 31 30282 
11. 00 184. 33 14. 76 20. 88 . 1055D+06 192. 17 26.10 
es 5 0 291.05 ess 23.04 » 8355D+05 105.79 de! 
200 282.49 17.08 60. 84 .4231D+05 13.42 Drs 
251) 0 oo) 229 TSea2 P3592 22o2.Ut05 = 205 1 ieeZ 
13. 00 36.98 2057 LEG 52 7309 7Dt05 = 6.506 15. 16 
a2 0 LOT 1 22:05 36. 00 . 8161D+05 7 On25 Sie oe 
14. 00 17 2e0 24.17 18. 36 7 13230506 9251.60 24.19 
14.50 250. 72 26. 56 22022 727100705 142.56 13.55 
15. 00 336. 19 Zone 105. 84 noo 20e05 =i. 77 13,36 
ib) 2 0 70. 18 Soe O 54.72 - 7805D+05 48.42 22709 
16. 00 173 63 35.48 15. 84 . 1645D+06 520.78 T9539 
16. 50 287.51 597.10 37.08 . 9863D+05 108. 50 11. 38 
00 D2. OS 43.07 59.04 - 21550405 59.16 222109 
i) 0 190.75 47.41 oo o2 . 1971D+06 705. 49 Ti 
18. 00 342. 34 52.14 69.48 .9377D+05 50. 26 20. 66 
18. 50 148. 81 oe, 12. 60 . 2192D+06 776. 24 9. 76 
19. 00 331.41 62782 53. 28 ~1295D+06 160. 46 16. 76 
9.50 171. 43 68. 81 10. 80 . 2601D+06 990. 18 7 oae 
20. 00 30°23 1320 46.44 . 1647D+06 29109 oir 
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D. BLADE ELEMENT RESULTS FOR BLADE ONE 


BLADE 
POSTION 
(DEG) 


54. 48 
56250 
66. 76 
79227 
95a? 
114. 25 
1350 
159. 43 
184. 89 
211236 
240. 10 
269. 43 
2 OL 2 
330. 89 
Zend 2 
55a 
69. 69 
104. 66 
140.91 
178. 69 
218. 28 
260. 03 
304. 34 
351. 66 
42.49 
97.40 
156. 98 
221791 
29290 
10272 
96.19 
196219 
Z293sG0 
47.51 
172386 
310m 
102. 34 
268. 81 
91.41 
291. 43 
150. 24 


147. 
140. 
or 
TY. 
75; 
46. 
a2 
ZS: 
ei. 
Zi 
36. 
hele 
140. 
154. 
154. 
15%: 
122. 
oe 
a2. 
25. 
25. 
54. 
144. 
154. 
147. 
64. 
Lae 
Zoe 
1 2c. 
154. 
ole 
Ze 
113: 
136. 
21. 
136. 


43. 
54. 
82. 
ZN; 


NDNDWONNVAWAAMWAANW WODWNHAWOAWOODOANEAVAENHVAWAMDMFNVWOAANHLOAAN HFA 


Z 


147. 
140. 
129. 
108. 
72s 
46. 
a2: 
ye 
2 
21. 
Se 
12: 
3 Gi. 
11. 
154. 
147. 
PLS: 
54. 
Zo: 
ois 
2D: 
46. 
140. 
154. 
144. 
Di 
2D: 
EDs 
ela 
11: 
DO. 
21. 
oS. 
22: 
18. 
A. 
ea 
nS 
50. 
43. 
18. 


ON ON KF DHDOFAKDENYVUA NN HADOWE ANN WOWAAAONWOFLFAANHLAOAAODOALA 


3 


147. 
140. 
129; 
108. 
2. 
43. 
ZGe 
21; 
206 
Ze 
326 
64. 
132: 
Lae 
154. 
147. 
11> 
ir 
28. 
Ze 
25: 
43. 
13 
15e 
140. 
50. 
213 
ale 
9: 
144. 
43. 
18. 
68. 
104. 
bar 
82. 
Zoe 
18. 
2>: 
aor 
14. 


FNDBDONDWOF FON OAAKAHLENYVWNHONNHNHWENVAODANNHWOHLFAAAMWMNHOOALA 


ANGLE OF ATTACK 


( DEGREES) 
SECTIONS 

& > 6 7 8 9 10 
147.6 147.6 147.6 144.0 144.0 144.0 144. 
140.4 140.4 140.4 140.4 136.8 136.8 136 
129.6 126.0 126.09126.0 122745122745 
104.4 100.8 100.8 97.2 93.6 93.6 90 
68.4 64.8 64.8 61.2 57.6 57.6 54 
43.2 39.6 39.6936.0 362500036, 00mee 
28.8 (26.8 25.2525. 2 (2575 eee 
21.6 2.6 Zis60e21.6 W820 cr eee 
21.6 18.0 18.0 5128.0 18.0 18. 000s 
18.0 18.0 18.0 18.0 18. 0951474 
28.8 28.8 25.2 025.2 2502552) 300g 
57.6 50.4 46.8 43.2 39.6 36,0522 
129.6 122.4 118.8 111.6 104.4 97.2 90 
151.2 147.6 144.0 144.0 140.4 136.8 133. 
151.2 151.2 147.6 147.6 147.6 144.0 140. 
144.0 144.0 140.4 136.8 133.2 129.6 126. 
108.0 1900.8 97.2 90.0 82.8 75.6 
46.8 43.2 39.6 39.6 36.0 36.70 gam 
28.8 25.2 25.2 °°25.2 213092) 60geee 
21.6 21.6 21.6 921.6 18.05 1820 Rare 
25.2 @ie6 21:6 921.6 1820 Sale 30 sacs 
39.6 36.0 32.4 926.8 25. 2.022552 fae 
126.0 118.8 108.0 97.2 86.4 57902 Wag 
151.2 147.6 144.0 140.4 136.8 129.6 126. 
133.2 129.6 122.4 115.2 104. 49897°2 Bees 
46.8 39.6 36.0 32.4 32.49925.68 Waa 
21.6 21.6 18.0 18.0 1850 "16.0 ike 
21.6 18.0 18.0 18.0 14.4 14.4 14. 
75.6 57.6 46.8 3976 32384 928.6 ee 
140.4 13322 122.4,11126 9752 82.5 |v 
36.0 32.4 28.8 72532 21269°21.6 fie 
18.0 14.4 14.4 14.4 14.4 10.8 10. 
46.8 36.0 28.8 25.2 21.6 1830 de 
82.8 64.8 50.4 39.6 32.4 28358 Ze 
144.4 14.4 14.4 10.8 10.8 10.8 10. 
54.0 39.6 28.8 21.56 18.0 98.0 de 
21.6 18.0 18.0 14.4 14.4 10.8 10. 
14.4 10.8 10.8 10.8 752 Bad 7 
21.6 18.0 14.4 14.4 10.8 10.8 10. 
16.0 14.4 10:8 516.8 10:8 Uae 7 
10-3 31056 ° 10a ted Tare die 2 7 


60 


WHOM OF ANF WOON FL FWOLOEFADOAHLODOFNHNOFAHLOONFLOCHLWO 


SECTION VERTICAL LIFT (LBS. ) 


POSITION 


BLADE 
(DEG) 
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i | anna ae | aoNre ra ~t NAN IMMOMM WO 
+ — rc eon! 
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t ot ac ae aN a NO rei Or 
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ANd! AKRANDS EMOTE NOTON ITH ANON “MIND © FOND 
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APPENDIX B. CH-46 ROTOR BLADE PROPERTIES 


Blade Station Pitch CG Flapwise Chordwise Torsional 
Station Radius Inertia Offset Stiffness Stiffness Stiffness 
(r/R) (in) (lb. *s*s) ein) Cibo in) (1b. inin) Clb: in*in) 
OeO16g7 5.5 Sa, 0.0 570. OE06 285. OE06 70. OE06 
0.05 4.85 oe 0.0 160. OE06 285. OE06 70. OE06 
0.10 6.90 nol 0.0 81. OE06 345. OE06 70. OE06 
0.15 erlS Lee 0.028 160.0EO6 #£260. OE06 67. OE06 
0. 20 0.78 lees 5 04 48. OE06 218. OE06 28. OE06 
O25 .65 2690 -.93 20. OE06 370. OE06 18. OE06 
0.30 soe Zell 69 21.0EO06 420. OE06 17. OE06 
0.35 ~505 2,09 74 20.0E06 412. 0E06 17. OE06 
0.40 ~ 499 2206 -.78 20.0EO06 405. 0E06 17. OE06 
0.45 - 48 2.04 -. 82 19. 5E06 395. OE06 16. 9E06 
0.50 ~45 2. 10 -. 84 19. 5E06 392. OE06 16. 5E06 
Oe 55 ~ 42 2, 00 -, 86 19. 5E06 390. OE06 16. 2E06 
0.60 - 42 Zale -. 58 16.0E06 410. 0OE06 15.5806 
0. 65 ~ 42 Zl -.51 16.0E06 420. OE06 15. 2E06 
0.70 ~ 42 2 lel -.51 16.0E06 420. 0E06 15. 2E06 
O75 ~ 42 2 -. 76 16.0EO06 570. OE06 15. 2E06 
0. 80 ~41 2 ale = 1 16.0E06 420. OE06 15. 2E06 
O.S5 ~45 2. 12 -, 34 18.0E06 435. 0E06 16. 5E06 
0.90 510 2.40 25 15.0E06 440. OE06 16. OE06 
0.95 oo 2555 33 12.0E06 420. 0E06 14. OE06 
1. 00 1. 30 3.010 -.63 18. OE06 755. OE06 15. 2E06 
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APPENDIX C. GENHEL FLOW DIAGRAM 
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StrAve 
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RO ton 
Ognere 


APPENDIX D. MATLAB META FILES 
A. U-COMPONENT MODEL FOR WIND TUNNEL 


Mmi=( oo.) =-.2 -3.2 -7.2 -8. 2] *10 
fu=[0 3/1000 35/1000 600/1000 999/1000 1] 
[ bu, au] =yulewalk(6, fu,mu) 

[ Au,Bu,Cu,Du] =TF2SS( bu, au); 

[ Zu, Polesu, Ku] =SS2ZP(Au, Bu,Cu,Du, 1) 
x=rand(1,4000); 

yu=filter(au,bu,x); 

pack 

Pu=spectrum(x,yu, 200,100); 
disp('Hit RETURN to see plots’); 
pause 

meta comp_u 

spcplotu(Pu,1000); 


B. V-COMPONENT FOR WIND TUNNEL MODEL 
clear 

mv=[0 .1 *-.7 -4.7 -7.2 -10. 2] *( 10/1. 73) 
fv=[0 3/1000 35/1000 600/1000 999/1000 1} 
[ bv, av] =yulewalk(6,fv,mv) 

[ Av,Bv,Cv,Dv] =TF2SS(bv, av); 

{ Zv,Polesv,Kv] =SS2ZP( Av, Bv,Cv,Dv,1) 
x=rand(1,4000); 

yv=filter(av,bv,x); 

pack 

Py=spectrum(x,yv,200, 100) = 

disp('Hit RETURN to see plots'); 

pause 

meta comp_v 

Bee loty( Pv, 1000); 


C. W-COMPONENT FOR WIND TUNNEL MODEL 


mw=([O .1 -.3 -2.7 -4.5 -6.5]*10 
fw=[0 3/1000 35/1000 600/1000 999/1000 1} 
[ bw, aw] =~yulewalk(6, fw,mw) 

[ Aw, Bw, Cw, Dw] =TF2SS( bw, aw); 

[ Zw, Polesw, Kw] =SS2ZP( Aw, Bw, Cw, Dw, 1) 
x=rand(1,4000); 

yw=filter(aw,bw,x); 

pack 

Pw=spectrum(x,yw,200, 100); 

disp ('Hit RETURN to see plots'); 
pause 

meta comp_w 

spcplotw( Pw, 1000); 
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D. U-COMPONENT FOR PROTOTYPE MODEL 


mu=[1 .1 *.2 -3.2 -7.2 -8. 2] *10; 
fu=[ 0 3/1000 35/1000 6000/1000 999/1000 1]; 
[ bu, au] =yulewalk(6, fu, mu) 

[ Au,Bu,Cu, Du] =TF2SS(bu, au); 

[ Zu, Polesu, Ku] =SS2ZP( Au, Bu, Cu, Du, 1) 
x=rand(1,4000); 

yu=filter(au,bu,x); 

pack 

Pu=spectrum(x, yu, 200, 100); 
disp('Hit RETURN to see plots’); 
pause 

meta comp_4u 

speplt4u( Pu, 1000/14); 


E. V-*COMPONENT FOR PROTOTYPE MODEL 


mv=[O0 .1 °.7 -4.7 -7.2 -10.2}*(10/1. 73); 
fv=[0 3/1000 35/1000 6000/1000 999/1000 1]; 
[ bv, av] =yulewalk(6, fv,mv) 

[ Av, Bv, Cv, Dv] =TF2SS( bv, av); 

[Zv, Polesv, Kv] =SS2ZP( Av, Bv, Cv, Dv,1) 
x=rand(1,4000); 

yv=filter(av,bv,x); 

pack 

Pv=spectrum(x,yv,200,100); 

disp('Hit RETURN to see plots’); 

pause 

meta comp_4v 

spceplt4v(Pv,1000/14); 


F. W-COMPONENT FOR PROTOTYPE MODEL 


mw=[0 .1 -.3 -2.7 -4.5 -6.5]*10; 
fw=[0 3/1000 35/1000 600/1000 999/1000 1]; 
[ bw, aw] =yulewalk(6, fw,mw) 

[ Aw, Bw, Cw, Dw] =TF2SS( bw, aw); 

[ Zw, Polesw ,Kw] =SS2ZP( Aw, Bw, Cw, Dw, 1) 
x=rand(1,4000); 

yw=filter(aw,bw,x); 

pack 

Pw=s pectrum(x,yw,200, 100); 

disp ('Hit RETURN to see plots’); 
pause 

meta comp_4w 

spcplt4w(Pw,1000/14); 


f is a vector of frequency points, specified in the range between 0 and 1, where 1.0 
corrsponds to half the sample frequency (the Nyquist frequency). 


m 1s a vecotr containing the desired magnitude response at the points specified in f 
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b and a are row vecotrs containing the n + I coefficients of the order n filter whose 
frequency-magnitude characteristics match those given in f and m. 


[A,B,C,D] are the state space vectors representing the digital filter transfer function. 


[Z,P,K] are the vectors representing the minimum phase zeros, poles and filter gain 
of the digital filter transfer function. 


BB 


APPENDIX E. INPUT FOR DYSCO MODULES 


H46A-1 /MODEL ON FILE U1 
Ke tetetetetedek keictetetetetekktetetek fete MODEL H46A-1 tektcktkkkkicteksekiekikdsedeteteteteiexk 
H46 BASELINE 
INDEX COMP NO. DATA SET FORCE DATA SET 
1 CRE3 2  H46A2 FRA3. WIND 
REQUIRED DS/DM= NACA0012/AIRFOIL 
2 CCEO 2 AFT_HEAD NONE 
3 CFM2 1 H46 NONE 
4, CSF1 GEAR FSS1 WAVES 
5 cho COUP_G NONE 
6 CES1 DROOPSTO NONE 
7 CSF1 BLADE_WG NONE 


FetedetetekeTeHRHRKRR A eK RK RRC KIRK KE ARK IK ERR AKR ERR RRREKIRRIERREKIKKKK KRIK ET 


HARRI etek tek KKK IKK RIKI IK KIRKE IK IKEA IK KTR REIKI RRRRRKEKRKKK I 


1 VSOUND 
2 RHO 


GLOBAL VARIABLES 


- SOUND VELOCITY = 1.11700E+03 


- AIR DENSITY RATIO = = 1. 00000E+00 
PAA A A HAA HAHEI KHER IER IRATE IIR IIA IRA ISIA IN IRI AA IR IIA ITAA RAI IN 
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Hake RK H46A2 /CRE3 aK RAK 


H46 AFT ROTOR (3 BLADES: B2) 


Vetere reset tee tek Keke Kt deke terete RK KR EKER KREUK RR 


INPUT FOR ROTOR COMPONENT CRE3. ROTOR ELASTIC BLADES 
1 JV - INPLANE DOF = YES 
2 JW - OUTPLANE DOF = YES 
3 JP - TORSION DOF = YES 
4 JS - SHAFT PERTURBED DOF = NO 
5 IX - XHUB( LONG) DOF = NO 
6 JY - YHUB(LAT) DOF = NO 
De az - ZHUB( AXIAL) DOF = YES 
8 JAX - ALFX(ROLL) DOF = NO 
9 JAY - ALFY(PTCH) DOF = YES 
10 JAZ - ALFZ(YAW) DOF = NO 
11 NV - NO. OF INPLANE MODES= 1 
12 NW - NO. OF OUTPLNE MODES= 2 
13 NP - NO. OF TORSION MODES= 1 
14 NB - NO. OF BLADES = 3 
15 PHL - PTCH HORN LNGTH (IN)= _ 8. 67000E+00 
16 PHSTA  -- PTCH HORN STA (IN) = 4.30000E+01 
17 NX - NO. OF STATIONS = 20 
18 ITYP - MODE INPUT 10R2 = 2 
19 X - (REAL) STATIONS 
5.11000E+00 1.42000E+01 4.30000E+01 5. 80000E+01 
7.65000E+01 9.18000E+01 1.07100E+02 1. 22400E+02 
1.37700E+02 1.53000E+02 1.68300E+02 1.83600E+02 
1.98900E+02 2.14200E+02 2.29500E+02 2.44800E+02 
2.60100E+02 2.75400E+02 2.90700E+02 3. 06000E+02 
20 NIP - INPLANE HINGE STA = 2 
21 NOP - OUTPLANE HINGE STA = 1 
22 NTOR - PTCH BEARING STA = 3 
penctPP - IP MODAL DAMPING = 0.00000E+00 
24 COPP - OP MODAL DAMPING = 0.O00000E+00 0. 00000E+00 
25 CTORR  - TORSION MODAL = 0.00000E+00 
26 IBIP <r BO ORs? = 1 
27 IBOP - OP BC 1 OR 2 = fl 
28 IBTO - TORSION BC 10R2 = 2 
29 NI - NO. OF IMPLICIT DOFS= 9 
30 NIDOF  - (DOF) IMPLICIT DOF NAMES 
NGOh NGI WEP “amever 12 WGr — 22 
WTIP 2 WGT 13 WGT 23 WTIP 3 
31 XSTA - (REAL) STAS FOR IMPLCT DOFS 
9.83200E+01 1.95600E+02 3.06000E+02 9.83200E+01 
1.95600E+02 3.06000E+02 9.83200E+01 1.95600E+02 
3. O6000E+02 
32 DIST - (REAL) DISTANCE TO EC 
0. 00000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 
0. OOOOOE+00 
33 IBN - BLADE NO. FOR IDOFS 
1 1 1 2 2 
2 3 3 3 


a2 


34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 


50 


51 


a2 


Bye 
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a0 


SE 


SEA 


KM1 


KM2 


URE 


—mo coqcmoiti oo°ocncoeo 1oc$c’tooo .  ece oc 


' 
>) 


we 
=2 


IP SPRING RATE = Q.00000E+00 
IP DAMPING RATE =" 7. 85380E+03 
OP SPRING RATE = Q.  O00000E+00 
OP DAMPING RATE = 0. 00000E+00 
TORSION SPRING RATE = 3.57793E+03 
TORSION DAMPING RATE= 0. Q0O000E+00 
RPM = 1.70770E+01 
ROTATION DIRECTION = =1 
AZIMUTH OF REF BLADE= 7.92710E+01 
PRECONE ANGLE (DEG) = 0. 00000E+00 
HUB WEIGHT (LB) = 0. QO0000E+00 
HUB M.O.I. ABOUT Y- = 0. 00000E+00 
ROOT PTCH ANG (DEG) = 3. 00000E+00 
NONLIN TERMS = NO 

UNIFORM BLADE = NO 


(REAL) MASS PER UNIT LENGTH 


. 50000E+00 5.5Q0000E+00 1.00000E+00 8. 0Q0000E-01 
.50000E-01 5.20000E-01 5.05000E-01 4.99000E-01 
.80000E-01 4.50000E-01 4.20000E-01 4. 20000E-01 
.20000E-01 4. 20000E-01 4. 20000E-01 4. 10000E-01 
-50000E-01 5.QQ000E-01 8.50000E-01 1. 30000E+00 
(REAL) CG OFFSET FROM EA 
.Q00000E+00 0.Q0000E+00 2.50000E-01 2. 20000E-01 
. 30000E-01 -6.90000E-01 -7.40000E-01 -7. 80000E-01 
. 20000E-01 -8.40000E-01 -8.60000E-01 -5. 80000E-01 
. LOOOOE-01 -5.10000E-01 -7.60000E-01 -5. 10000E-01 
-40000E-01 2.S50000E-01 3.30000E-01 -6. 30000E-01 
(REAL) AREA CENTROID OFFSET 
- QOO00E+00 0. 00000E+00 0.00000E+00 0. 00000E+00 
. QOOOOE+00 0.00000E+00 0. 00000E+00 QO. 00000E+00 
. QOOO0E+00 0.00000E+00 0.Q0000E+00 0. 00000E+00 
-QOO000E+00 0.Q0000E+00 0.00000E+00 QO. 00000E+00 
-QO0O000E+00 0.OQ0000E+00 0. 00000E+00 0. O0000E+00 
(REAL) MASS ROG ABOUT 
-QQOO000E+00 0. 00000E+00 Q.00000E+00 0. 00000E+00 
-QOOO00E+00 0. 00000E+00 0.Q00000E+00 0. O0000E+00 
- OOO00E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
-QOO00E+00 0.00000E+00 0.Q0000E+00 0. 00000E+00 
-QOOO00E+00 0. 00000E+00 0.00000E+00 0. 00000E+00 
(REAL) MASS ROG ABOUT 
. OOOOOE+00 0. Q0000E+00 0.Q0000E+00 0. 00000E+00 
- QOOOOE+00 0. 00000E+00 0.00000E+00 0. 00000E+00 
- QOOO00E+00 0.00000E+00 0.0Q0000E+00 0. 00000E+00 
- QOOOQ00E+00 Q.00000E+00 Q0.Q00000E+00 0. 0Q0000E+00 
- OOOO0E+00 0.00000E+00 0.Q0000E+00 0. 00000E+00 
(REAL) AREA ROG OF CROSS 
. QOOO0E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
.QQOQ000E+00 0. Q00000E+00 Q.Q0000E+00 0. 00000E+00 
- QOOOO0E+00 0. 00000E+00 0.00000E+00 0. 00000E+00 
- OOOOOE+00 0.Q0000E+00 0. Q0000E+00 0. 00000E+00 
-QO000E+00 0. O00000E+00 0.00000E+00 0. 00000E+00 
(REAL) PRETWIST RATE DEG/IN 
-QOO000E+00 Q.Q00000E+00 0.Q0000E+00 2.90200E-01 
- 80000E-02 -2. 80000E-02 -2. 80000E-02 -2. 80000E-02 
- 80000E-02 -2. 80000E-02 -2. 80000E-02 -2. 80000E-02 
- 80000E-02 -2.80000E-02 -2. 80000E-02 -2. 80000E-02 
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-2. 80000E-02 -2.80000E-02 -2.80000E-02 -2. 80000E-02 
Sey OG - (REAL) CHORDWISE EI*10E-6 
2.85000E+02 2.85000E+02 8.00000E+01 2. 75000E+02 
3. 7OOO0E+02 4. 20000Et+02 4.12000E+02 4. 05000E+02 
3.95000E+02 3.92000E+02 3.90000E+02 4. 10000E+02 
4.20000E+02 4.20000E+02 5.70000E+02 4. 20000E+02 
4.35000E+02 4.40000E+02 4. 20000E+02 7.55000E+02 
SI - (REAL) BEAMWISE EI*10E-6 
5. 7VOOOOE+02 1.60000E+O2 2.10000E+02 7. 00000E+01 
2,00000E+Omeee. LOOCUOETOL “200000RTOl 2. 00000E+01 
J OS000RtOmr!. SSODOEKTUL “1.95000 -OL >) lso0000E+0] 
1.60000E+01 1.60000E+01 1.60000E+01 1. 60000E+01 
IeSG0C0R+O hee. SOCOOETOT 1. 2Z20000E+TOl Tes0000E+01 
58 EA - (REAL) SECTION EA*10E-6 
0.O00000E+00 0.0Q0000E+00 Q0.O00000E+00 0. O00000E+00 
0. OOO00E+00 0.O00000E+00 0.O00000E+00 OQ. O00000E+00 
0.00000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0.OO000E+00 0.00000E+00 0.O00000E+00 0. 00000E+00 
0.00000E+00 0.O00000E+00 0.00000E+00 0. 00000E+00 
be GJ - (REAL) SECTION GJ*10E-6 
7.QQOO00E+0O1 7.Q0000E+01 7.00000E+01 3. 40000E+01 
PeOCOORTOIea. 7O0COE+TOL 1. 70000E+01 1. 70000E+01 
imo O00 +O. ©SO00E+01 1.6Z000E+O1” 1. 55000E+01 
oe OCCURS Otel) 2OUOR+O1. 1. 52000E+TOl- 1. 5Z000E+01 
I OS000E+O01 1, 60000ETO1 1.40000EtO1 1.52000E+01 
60 EB1 - (REAL) CROSS SEC INTEGRAL 
0. OOO00E+00 O.Q0000E+00 0. O00000E+00 0. O00000E+00 
0. OOOOOE+00 O.O00000E+00 Q.Q0Q0000E+00 QO. Q00000E+00 
0.00000E+00 0.00000E+00 0.O00000E+00 0. O00000E+00 
0. OOO00E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0.Q00000E+00 0.00000E+00 0.O00000E+00 0. 00000E+00 
61 EB2 - (REAL) CROSS SEC INTEGRAL 
0.Q0000E+00 0.00000E+00 0. 00000E+00 0. 00000E+00 
0. OO0000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0. OOOOOE+00 0.00000E+00 0.O00000E+00 0. 00000E+00 
0. Q0000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0. OOOOOE+00 0.00000E+00 0.00000E+00 0. 00000E+00 
62 EC1 - (REAL) CROSS SEC INTEGRAL 
0. OO000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0. 00000E+00 0.00000E+00 0.00000E+00 0. O00000E+00 
0. OOOOOE+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0.Q0000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0. OOOO0E+00 0.O00000E+00 0.00000E+00 0. 00000E+00 
63 EC1STA - (REAL) CROSS SEC INTEGRAL 
0.QO0000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0. Q0000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0. QOOO000E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0. QOOOOE+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0. OOOO00E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
64 JIL - INTERNAL LOADS = YES 
65 NXIL - NO. OF STATIONS = 20 
66 JIPIL - INPLANE MOMENTS = NO 
67 JOPIL - OUTPLANE MOMENTS = YES 
68 JTORIL - TWIST MOMENTS = 


FEMME He FoF Te Tee TE TC TE TC He HEHE TE He FETE eA He Fe FEN Fe FETE Te HTC ETE TE FETE TC TCE TE FETE TE TE TOPE TET TE TCT TTI TE FCF T CTI TC ICI ICI 
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Bait a aa ek AI H46 /CFM2 Fewer Wierd Hes 


H46 FUSELAGE 


FereFesewe sede se Fess He He WTF TH HN HH HIE KEK RRR EKER RRR EERE RRR KK KR TI 


INPUT FOR STRUCTURAL COMPONENT CFM2. MODAL FUSELAGE 
RBM - RIGID BODY MODES YES 
IXCG - LONGITUDINAL NO 
IYCG - LATERAL NO 
IZCG ~ VERTICAL YES 
IROLL - ROLL NO 
ERTCH =" PLICH YES 
TYAW - YAW NO 
CG CG STATION (IN) 2. 94500E+02 


Oe ee ee 
OVW Wr DAMN FWNHOKP OW WIA NHK WH He 


NMODE - NO. OF ELASTIC MODES 0 
NR - NO. OF ROTORS Z 
NROT ~ ROTOR NUMBERS 1 2 
XROT - ROTOR STATIONS 7.40000E+01 4. 76000E+02 
ZROT - ROTOR VERTICAL HT 8. OOOOOE+01 1. 28000E+02 
ASF - FWD SHAFT ANGLE 9. 50000E+00 7.00000E+00 
ASL - LAT SHAFT ANGLE 0. OOOOOE+00 0. 00000E+00 
IX - HUB TRAN DOF - LONG NO NO 
ley. - HUB TRAN DOF - LAT NO NO 
IZ - HUB TRAN DOF - AXIAL= YES YES 
TAX - HUB ANGL DOF - ROLL = NO NO 
IAY - HUB ANGL DOF - PITCH= YES YES 
21 IAZ ~ HUB ANGL DOF - YAW = NO NO 
22 NI - NO. OTHER IMPLCT DOF= 0 
23 MASSL - FUSELAGE MASS (LB) = 1.94400E+04 
24 IMYF - PITCH MOI ABOUT CG = 1.16700E+05 


FereaeTeseve dese ke Tere Hk Fe Were TWH ee FeFe WF HH HAM Te WH He FETCH Fe Fe TE HT IK FEM Te WIC TE PRT HT TRIER I I 
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deacdeiek reasekledcdetek AFT_HEAD/CCEO Haddad Tek ees 


CONTROL RODS FOR H-46 AFT HEAD 


Feeder seve vese Te seae Te severe Weve de Meese HeMe de He PeA AI MI MAA AIT AT AIH ATE HACIA AI TAI TAI TAIT CII I 
INPUT FOR CONTROL SYSTEM COMPONENT CCEO. CONTROL RODS 


1 KROD - CONTRL ROD STIFFNESS= 3. 00000E+03 
Kicacacacadeaaksedesteskdesdedseacaede sede teak dese aksdevetededededevede ded dete devote deeded Heer PeeKRRKKRK KK KR 


Jennie deeded GEAR /CSF1 Ke I 


LANDING GEAR 


HHAA TAIN IRIE RIK IK TRIER ITER TEEN IERIE IER TEETER IER IER IERIE IKEA IK IKARIA ARIK 
INPUT FOR COMPONENT CSF1. FINITE ELEMENT 


1 NCDF - NUMBER OF DOF = eZ 
2 CDFLI - DOF NAMES = FG 0 AG 0 
Bactl - (REAL) MASS MATRIX 
NULL MATRIX 
4 CC - (REAL) DAMPING MATRIX 
NULL MATRIX 
5 CK - (REAL) STIFFNESS MATRIX 


DIAGONAL MATRIX (DIAGONAL VALUES PRINTED) 


6.50000E+03 6. 50000E+03 
6 CF - FORCE VECTOR = Q0.00000E+00 0. 00000E+00 


7 IGR - GLOBAL REFERENCE 
HAAR KATE KAIRIE IRENE TICK HICH IK TINE AIK TRAE RI IK IKI IIR IKARIA 
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Hered deka aaa COUPSG ay Olen PICKIN WII IT 


RELATE GEAR DOF TO FUSELAGE 


HIKAANAAKAKAK RIKI K NINN IN KIKI KIRIN INI IKK IN IK IRIAN INI IN I IAI 
INPUT FOR COMPONENT CLC1. LINEAR CONSTRAINTS 


1 NCDF - NUMBER OF DOF = 2 
2 CDFLI - DOF NAMES = ZCG 1000 PTCH1000 
3 NCIDF - # OF CONSTRAINT EQNS= 2 
4 CIDFLI - IMPLICIT DOF NAMES = FG 0 AG 0 
> COEF - (REAL) COEFFICIENT MATRIX 
GENERAL MATRIX 
ROW 1 
1. OOOOOE+00 2. 64000E+02 
ROW 2 


1. OOOOOE+00 -5. OOOOOE+01 
HIKAAKKA KA KAAIKIN KIER II IN IAAI IH IAAI AIK AIR I NICK TIA IH IAA FAIA HC II 


DROOPSTO/CES1 ON FILE U1 


tok terektedete dae tease DROOPSTO/CES1 nists deteaciek tik 


DROOP STOP FOR ROTOR HEAD (FLAP STOP) 


KikkkakdkktkkickiktkkkackkekkckkkKkKeKkkckkKKkckkekkeKKeK KKK KKK RK KKK KK KK 
INPUT FOR COMPONENT CES1. ELASTIC STOP 


1 MCDF - # OF DOF-EXCEPT BASE= 6 
2 CDFLI - (DOF) DOF NAMES 
OP ZITO OP 2120 OP «2210 OR 3z220 WOR 
OP 2320 
3 BASE - EXISTNCE OF BASE DOF= YES 
4 CDFLBI - BASE DOF NAME = ZHUB 1000 
561 - UPPER DAMPING COEFF = Q.00000E+00 
6 C2 - LOWER DAMPING COEFF = 0. 00000E+00 
ao - UPPER SPRING COEFF = 0.00000E+00 
6 KZ - LOWER SPRING COEFF = 0. 00000E+00 
SODEDT - UPPER GAP SIZE = 2.10000E-02 
10 DELT2 - LOWER GAP SIZE = 0.00000E+00 


Kick KIcterk HI Keke WICH He Fe Fe Fe FEK MAE FeFCT FIER ACH RAEI FETE HK A FE FER PKA RIK KIKK KIKI Fe KIRK IK 
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Hitter ee BLADE _WG/CSF1 KeK KKK IK 


Blade weight DOF for CRE3 


HHA HI KT TIT ITT TRI RTI IEEE IRR RRR TERRIER 
INPUT FOR COMPONENT CSF1. FINITE ELEMENT 


1 NCDF 
2 CDFLI 


3 CM 


4 CC 


> CK 


6 CF 


ROW 
ROW 
ROW 
ROW 
ROW 
ROW 
ROW 
ROW 
ROW 


- NUMBER OF DOF = 9 

- (DOF) DOF NAMES 
WGT 1 WGT 2 WIIP O WGT 12 WGT 22 
Withee e2 WClee 15 7WGl o23eeWwIT Pe. 3 

- (REAL) MASS MATRIX 
NULL MATRIX 


- (REAL) DAMPING MATRIX 
NULL MATRIX 


- (REAL) STIFFNESS MATRIX 
NULL MATRIX 


- (REAL) FORCE VECTOR 
0. OOO0OO00E+00 0.00000E+00 0.00000E+00 0. 00000E+00 
0.O0000E+OO0 0. O00000E+00 Q.O00000E+00 0. O0000E+00 
0. OOOOOE+00 
- GLOBAL REFERENCE = YES 
- (REAL) LOCAL DOF VECTORS 

GENERAL MATRIX 


1 


0. OOOOOE+00 0. Q0000E+00 1. 00000E+00 


RO 


Q. QOOO00E+00 0.00000E+00 1. 00000E+00 


Ww 


0. 00000E+00 0.00000E+00 1. 00000E+00 


> 


0.O0000E+00 0.00000E+00 1. 00000E+00 


o,) 


0. OOOOOE+00 Q.Q00000E+00 1. 00000E+00 


ON 


0. OOOOOE+00 0. 00000E+00 1. 00000E+00 


7 


0. OOOOOE+00 0.00000E+00 1. 00000E+00 


8 


0. OOOOOE+00 0.00000E+00 1. 00000E+00 


9 


0. OOOOOE+00 0. 00000E+00 = 1. 00000E+00 


Fetetete RIK Teter te Ke Kee Ke Fete FeAe HK FFI TC TCR TCR Te PRTC KICK Te FCT ICTRACKER KET TREREETRIK TT RICK Fe 


8] 


Rake KKK 


WIND /FRA3 
LOW REYNOLD'S NUMBER AERODYNAMICS FOR NACA0012 AIRFOIL 


REQUIRES DS/DM NACA0012/AIRFOIL 
NO SEQUENTIAL FILES REQUIRED 


HIRAI HIT IAIN INN ITI IIIT IN IIIA KAIRIE I 
INPUT FOR FORCE FRA3. ROTOR AERO GENERAL 


TEQS AERODYNAMICS BY EQS = NO 
INFTAB INDUCED VEL BY TABLE= NO 
IUNSTD UNSTEADY AERO = NO/YES 


(REAL) WIND VELOCITY 


oe 3emaoaoas 


1 
a 
5 
4 VAIRH 
5 
6 
7 
8 


.22335E+01 0.00000E+00 2.47238E+01 
ASTALL STALL ANGLE (DEG) = 1. 20000E+01 
RFCT INDUCED VEL FACTOR = _ 1. 65000E+00 
TIPLOC TIP LOSS COEFFICIENT= 9.50000E-01 
XH - HUB EXTENT (IN) = 1.16000E+01 
9 ALT - VEHICLE HEIGHT (FT) = 0. 00000E+00 
10 K27 - TIP VORTEX COEFF = QO. O00000E+00 
11 CDO - BLADE DRAG COEFFAT = 1. 00000E-03 
TZ 5OtG - Q1C COEFFICIENT = 1. 00000E+00 
132-QZC - Q2C COEFFICIENT = 5. Q00000E-01 
14 ALAMDA - NONDIM INDUCED VEL = _ 0. 00000E+00 
15 NXA - NO. OF STATIONS = 10 
16 XAERO - (REAL) NONDIM AERO STATIONS 
1.Q00000E-01 2.00000E-01 3.00000E-01 4. Q0000E-01 
5.Q00000E-01 6.0Q0000E-01 7.Q0000E-01 8.0QQ000E-01 
9. QOOOO00E-01 1. 00000E+00 
17 NUMAF - NO. AIRFOIL TABLES = i 
18 AFTAB1 - NAME AF TABLE 1 = NACA0012/AIRFOIL 
19 NUMAF1 - NO. OF STATIONS AF 1= 10 
20 STA-AF1 - STATIONS FOR AF 1 
1 2 3 a 5 
6 7 8 9 10 
21 XACC - (REAL) A/C OFFSET FROM E.C. 
0. Q0000E+00 Q0.O00000E+00 0.Q00000E+00 0. 00000E+00 
0. OOOOOE+00 0.00000E+00 0. O00000E+00 0. 00000E+00 
0. OOOOOE+00 0. O0000E+00 
22 CHORDC - (REAL) CHORD (IN) 
0. OOOOOE+00 1.10000E+01 1.87500E+01 1.87500E+01 
1.87500E+O1 1.87500E+0O1 1. 87500R4FOilee, 17675008701 
1.87500E+01 1.87500E+01 
23 NX - NO. AERO FACTOR STAS= 10 
24 XF - (REAL) NONDIM FACTOR STAS 
1.Q00000E-01 2.00000E-01 3.00000E-01 4. Q0Q000E-01 
5. Q00000E-01 6.00000E-01 7.Q0000E-01 8. 0O0Q00E-01 
9. OOOO0E-01 1. 00000E+00 
25 FL - (REAL) FACTORS FOR CL 
1. QO000E-01 5.00000E-01 7.00000E-01 1. 00000E+00 
1. OOOOOE+00 1.00000E+00 1.00000E+00 1. 00000E+00 
9.80000E-01 8.90000E-01 
? 26 FD - (REAL) FACTORS FOR CD 
1.Q00000E-01 5.00000E-01 7.00000E-01 1. 00000E+00 
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-QOOO00E+00 1.00000E+00 1.Q00000E+00 1. 00000E+00 
. OOOOOE+00 1. 00000E+00 
(REAL) FACTORS FOR CM 
-QOOO00E+00 1. 00000E+00 1.00000E+00 _ 1. 00000E+00 
. OOOOOE+00 1.00000E+00 1.00000E+00 1. 00000E+00 


1. OOOOOE+00 1. 00000E+00 
FRAAKAK A IIHR INNATE TRIKE RIKEN IIE THINK IKE RIN ICTR RI IK IIR IIR IK RIE IKI HI I 


Za FM 


es a 
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APPENDIX F. VELOCITY FIELDS OF THE ROTOR BLADES 


(RHOADES) 


ALL VELOCITIES IN M/S 


o° VR 
RWS UY SKEW V 


2B nam 
Bog oge 
288 GBH 
Bont es 
yboxe 

Eécd%s 
s2e28E 
posers 
BR eee 
Bon iey 


0. 164 


-2.7356 
-3.438 
-3.766 
-3.055 


0.047 
oO. 


0.088 
0.024 
0.037 


0.1735 
-0.001 
-0.003 


-0.005 
-0.031 


HEARN V 


ete 
TERE 
pecce® 
zaRG oh 
Regare 
=Be0 

© ooo 
~SRree 
Eacedd 
BRT NOF 
: 


30° YR 


0.335 
38.584 
-3.859 

2.187 


0.449 


TURB W 


~0.64 
-3.514 
0.017 
0.084 


SKEW W 
-0.076 


-1.08 
4.2395 
-20.756 


-1.631 


TURS V 


-2.871 
—4.371 
0.034 
0.261 


SKEW V 
-0.016 


0.0937 
0. 

0.051 
0.237 
0.163 


Ris ¥ 


-0.002 
-0. 158 
-0.034 


MEAN Vv 
-0.069 
-0.025 


PET OSH 

BRe5 
Boccc® 
= OSE = 
Resyec 
-e8hq- 
Boe ad? 
-oRER 
enoo = 
a oe 
: 


ZO UPRS 
ANN RH 
Bonyoe 
ag2ahe 
Gocade 
zeroe 


“1.51 
1.458 


TURB V 
~0.917 
i. 
-0.225 —4.893 


0.16 
0.631 
0. 


SKEW V 
-0.015 


BRREL 


po 
geome 
: 


=PEELE 
ND OM 
BSaaee 
azarnnwet 
geeoks 
ooooo 
neR aR 
Ecosse 
2enVin @ 
gRBAKS 
yoogoe 
ee eae 
GST ade 
-ERSES 
bodice 


=SanY, 
keggte 
DEB KSR 
Focddc 
>ENORAR 
punene 
aie eS 
Byer 
c 


poset 


cues 


geet: 


>PonEes 


ON EOE S 
- = 


Fe scce 


pNoaes 


wn = 


rename 
¥ 
< 


POST T1046 


jae 
aNebon 
fA 


“O99 


-PEBee 


gocdc 


84 


APPENDIX G. LOOKUP TABLES FOR DYSCO FRA3 MODULE AND 
BLADE ELEMENT PROGRAM 


MACH NO. AOA LIFT DRAG § MOMENT 
(DEG) COEFF. COEFF. COEFF. 
0.060  -180.00 0.000000 0.022000 0.000000 
0.060 -176.40 -0.091100 0.049100 0.315800 
0.060 -172.80 -0.035400 0.087300 0.417300 
0.060 -169.20 0.152700 0.148200 0.384600 
0.060 -165.60 0.367900 0.231100 0.305300 
0.060 -162.00 0.570700 0.266200 0.290900 
0.060 -158.40 0.734400 0.353300 0.315100 
0.060 -154.80 0.863200 0.475100 0.349400 
0.060 151.20 0.952600 0.604700 0.383200 
0.060 -147.60 0.997300 0.740700 0.415400 
0.060 -144.00 1.013100 0.881500 0.445000 
0.060 -140.40 1.031500 1.025800 0.470800 
0.060 -136.80 1.043500 1.166400 0.491700 
0.060 -133.20 1.039100 1.293900 0.507100 
0.060 -129.60 1.008400 1.408400 0.519200 
0.060 -126.00 0.951200 1.510200 0.528500 
0.060 -122.40 0.881500 1.599400 0.535100 
0.060 -118.80 0.801000 1.676900 0.539200 
0.060 -115.20 0.710900 1.751000 0.541000 
0.060 -111.60 0.612300 1.822000 0.540600 
0.060 -108.00 0.506600 1.888200 0.538100 
0.060 -104.40 0.394800 1.947900 0.533700 
0.060 -100.80 0.278400 1.994300 0.527600 
0. 060 “97.20 0.158400 2.020500 0.519800 
0. 060 -93.60 0.036200 2.027700 0.510600 
0. 060 -90.00 -0.087100 2.022000 0.500000 
0. 060 -86.40 -0.210200 2.007900 0.488400 
0. 060 -82.80 -0.332000 1.985500 0.476000 
0. 060 -79.20 -0.451100 1.954300 0.462800 
0. 060 -75.60 -0.566400 1.915700 0.448700 
0. 060 -72.00 -0.676600 1.870300 0.433600 
0. 060 -68.40 -0.780500 1.817600 0.417500 
0. 060 -64.80 -0.876900 1.756600 0.400400 
0. 060 64.80 0.887100 1.757600 -0. 400400 
0. 060 68.40 0.795200 1.818000 -0. 417500 
0. 060 72.00 0.696100 1.870100 -0.433600 
0. 060 75.60 0.591000 1.915500 -0. 448700 
0. 060 79.20 0.480800 1.954300 -0. 462800 
0. 060 82.80 0.366600 1.985500 -0. 476000 
0. 060 86.40 0.249400 2.007900 -0. 488400 
0. 060 90.00 0.130300 2.022000 -0. 500000 
0. 060 93.60 0.010200 2.027700 -0.510600 
0. 060 97.20 -0.109700 2.020500 -0. 519800 
0. 060 100.80 -0.228400 1.994300 -0.527600 
0. 060 104.40 -0.344900 1.947900 -0. 533700 
0.060 108.00 -0.458100 1.888200 -0. 538100 
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APPENDIX H. WAKE TURBULENCE MODELING FILTERS FOR 330° 
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